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ABSTRACT 

We present an X-ray and multiwavelength study of 33 weak emission-line quasars (WLQs) and 18 quasars 
that are analogs of the extreme WLQ, PHL 1811, at z r; 0.5-2.9. New Chandra 1.5-9.5 ks exploratory obser¬ 
vations were obtained for 32 objects while the others have archival X-ray observations. Significant fractions of 
these luminous type 1 quasars are distinctly X-ray weak compared to typical quasars, including 16 (48%) of the 
WLQs and 17 (94%) of the PHL 1811 analogs with average X-ray weakness factors of 17 and 39, respectively. 

We measure a relatively hard (L = L16 ;!;q 32) effective power-law photon index for a stack of the X-ray weak 
subsample, suggesting X-ray absorption, and spectral analysis of one PHL 1811 analog, J1521H-5202, also in¬ 
dicates signihcant intrinsic X-ray absorption. We compare composite SDSS spectra for the X-ray weak and X- 
ray normal populations and find several optical-UV tracers of X-ray weakness; e.g., Ee II rest-frame equivalent 
width and relative color. We describe how orientation effects under our previously proposed “shielding-gas” 
scenario can likely unify the X-ray weak and X-ray normal populations. We suggest that the shielding gas may 
naturally be understood as a geometrically thick inner accretion disk that shields the broad line region from the 
ionizing continuum. If WLQs and PHL 1811 analogs have very high Eddington ratios, the inner disk could 
be significantly puffed up (e.g., a slim disk). Shielding of the broad emission-line region by a geometrically 
thick disk may have a significant role in setting the broad distributions of C IV rest-frame equivalent width and 
blueshift for quasars more generally. 

Subject headings: accretion, accretion disks - galaxies: active - galaxies: nuclei - quasars; emission lines - 
X-rays; galaxies 


1. INTRODUCTION 

Luminous X-ray emission is considered a universal prop¬ 
erty of active galactic nuclei (AGNs), and built upon this idea 
are extragalactic X-ray surve ys for finding AGNs efficie ntly 
throughout the universe (e.g.. lBrandt & Alexandedl2015L and 
references therein). Eor AGNs that are not radio loud (with 
a jet-linked X-ray enhancement) or X-ray absorbed, the 
X-ray-to-optical power-law slope parameter (aoxl3 has a 
highly significant correlation with 2500 A monochromatic lu¬ 
minosity (L^ ^nn A) across sa 5 orders of magnitude in UV lumi- 
nosity (e.g.. ISteffen et al.ll2006t llust et alJl2007l : iLusso et alJ 
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’2 oox is defined as oqx = “0.3838 ^//2 kev), with ^ and 
fi keV being the rest-frame 2500 A and 2 keV flux densities. 


120101) . X-ray emission from AGNs is believed to originate 
from the accretion-disk “corona” via Comptonization of disk 
optical/UV/EUV photons (e.g., [Turner & Milled 120091 and 
references therein), although the details of this mechanism re¬ 
main mysterious. 

Lew AGNs are found to be intrinsically X-ray weak, i.e., 
producing much less X-ray emission than exp ected from the 
aox-T 25 ooA relation (e.g.. iGibson et al.l 120081) . A few can¬ 
didates have been suggested recentl y based on Nuclear Sp ec- 
troscopic Telescope Array (AM5T4R; lHarrison et al.l2013l) ob¬ 
servations of sig nihcantly X-ray weak broad absorption l ine 
(BAL) quasars (iLuo et al.l 12013L I2014t iTeng et aT 120141) FI 
The best-studied intrinsically X-ray weak AGN is the type 1 
quasar PHL 1811 , a very bright {B = 13 .9) radio-quiet quasar 
at z = 0.19 (e.g., iLeighlv et alj|2007allbl) . It is X-ray weak 
by a factor of « 30-100 relative to expectations from the 
q;ox-T 25 oo a relation, and its X-ray weakness is likely in¬ 
trinsic instead of being due to absorption given its canon¬ 
ical X-ray spectrum (power-law photon index L = 2.3 ± 
0.1), lack of detectable photoelec tric X-ray absorption , and 
short timescale X-ray variability (ILeighlv et alJl2007bh . In- 
terestingly, PHL 181 1 also has an unusual UV spectrum 
(ILeighlv et alJ l20()7ah . which is dominated by strong Ee II 
and Ee III emission with very weak high-ionization lines. 
The rest-frame equivalent width (REW) of the C IV A1549 
line (6.6 A) is a factor of « 5 times smaller than that mea¬ 
sured from quasar composite spectra (30 A); this line is also 

BAL quasars ai'e identified by th eir broad (> 2000 Itm s“' wide) 
blueshifted UV absorption lines (e.g., IWevmann et alj 119911) : e.g., the 
C IV A1549 line. BAL quasars are in general X-ray wealc, often due to ab¬ 
sorption, but intrinsic X-ray wealtness is a viable explanation for a subset of 
BAL quasars. 
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blueshifted (by 1400 km s“') and asymmetric, indicative of 
a wi nd component in the broad emission-line region (BELR; 
e.g. . [Rich ards et a l]l201 1 ). Based on photoionization model¬ 
ing, suggested that many of the unusual 

emission-line properties of PHL 1811 are a result of the soft 
optical-to-X-ray ionizing continuum caused by the intrinsic 
X-ray weakness. 

The intriguing possible connection between the extreme 
emission-line and X-ray properties of PHL 1811 prompted a 
search for more such X-ray weak quas ars using UV emission¬ 
line selection criteria (iWu et al.ll201 IL hereafter Wll). A pi¬ 
lot sample of eight type 1, radio-quiet, non-BAL quasars with 
PHL 1811-like emission-line properties (including small C IV 
RLWs, large C IV blueshifts, and strong Le II and Le III emis¬ 
sion), termed PHL 1811 analogs, were selected for X-ray 
study. All of them turned out to be X-ray weak, by factors 
of > 4.8 to > 34.5, confirming the empirical link between 
the X-ray weakness and unusual UV emission-line proper¬ 
ties. However, an X-ray stacking analysis revealed a hard 
spectrum on average for this sample, albeit with a large un¬ 
certainty, suggesting that unlike PHL 1811 itself that appears 
to be intrinsically X-ray weak, these PHL 1811 analogs may 
often be X-ray absorbed. More PHL 1811 analogs selected in 
addition to the eight Wll pilot objects are clearly required to 
constrain better the nature of these extreme quasars. 

There is another small population of type 1 quasars 
that have X-ray and UV emission-line properties over¬ 
lapping with those of the PHL 1811 analogs; radio¬ 
quiet weak emission-line quasars (WLQs). Strong broad 
emission lines in the optical and UV are a characteris¬ 
tic featur e of radio-quiet quasar sIt was thus surpris¬ 
ing when iMcDowell et al.l (I l995h discovered the first WLQ 
PG 1407H-265, with unusually weak Lya, C IV A1549, 
C III] A1909, and H/3 lines. With the large spectroscopic 
quasar s ample provided by the Sloan Digital Sky Survey 
(SDSS; I York et al.l [200i, more WLQs were discovered. 
These were originally a t z > 2.2 where Lyg coverage is avail- 


able (r; 90 WLQs; e.s.jLan et al.lll999l;lAnderson et al.1 

o 

o 

iPlotkin et al.ll2008l UOlObl iDiamond-Stanic et al.ll2009 

), and 


later extended to l o wer redshifts (ft i 100 WLQs; e.g. 


Collin 

p et al.l 120051; IHrvniewicz et al.l 120101; IPlotkin et al. 

2010^ 

hi iNikolaiuk & Walterl 120121 iMeusinger & Balafkan 

2014|) requiring weak C IV and/or other lines at longer wave- 


lengths. The fraction of X-ray weak quasars among either 
the high-redshift or l ower-redshift W LQs is high (Ri 50%; 
IShemmer et al.|[200^ IWu et al.ll2012i hereafter W12), again 
suggesting a link between the weak UV line emission and X- 
ray weakness. 

Based on the overall similarities between the PHL 1811 
analogs and X-ray weak WLQs, Wll argued that PHL 1811 
analogs are a subset of WLQs, despite small technical dif¬ 
ferences in their UV line REW selection criteria^ WLQs 
contain both X-ray normal and X-ray weak quasars (we adopt 

In radio-loud systems, the line emission can sometimes be diluted by 
the synchrotron emission from a relativistic jet, as typically seen in BL Lac 
objects. 

PHL 1811 analogs were required t o have C IV REW < 10 A, while the 
WLQs in the Wll study were from the lPlotkin et all <2010bl) catalog which 
requires REW < 5 A for all UV emission features (e.g., C IV, C III], and 
Mg II). No apparent difference was found between the PHL 1811 analogs 
having < 5 A C IV REWs and those having 5-10 A C IV REWs. Therefore, 
the different REW criteria were considered a technical selection effect and the 
WLQ criterion could be relaxed to C IV REW < 10 A, which is the lower 3cr 
limit of the log-normal C IV REW distribution IDiamond-Stanic et alj|2009l : 
W12). 


an X-ray weakness factor of 3.3 as the dividing threshold be¬ 
tween X-ray weak and X-ray normal quasars; see Section|4T| 
below), while PHL 1811 analogs are likely X-ray weak WLQs 
due to the additional selection criteria of strong UV Le emis¬ 
sion and large C IV blueshift (W11; W12). A larger sample of 
PHL 1811 analogs than the pilot sample of eight would help 
examine further the above suggested connection. 

Except for the unusual UV emission-line and X-ray prop¬ 
erties, the PHL 1811 analogs and WLQs appear to be typical 
quasars in t erms of other obs ervable multiwavelength prop¬ 
erties (e.g., iLane et akl 120111 ; Wll; W12). Various expla¬ 
nations have been proposed for the nature of the PHL 1811 
analogs and WLQs, such as an anemic BELR where there is 
a significant deficit o f line-emitting gas in the BELR (e.g., 
IShemmer et akll^lOl) . a brief evol utionary stage where th e 
BELR is not fully developed (e.g.. iHrvniewicz et al.ll201(ih . 
a soft ionizing spectral energy distribution (SED) produced 
by the cold accretio n disk of a yery massiye black hole (e.g., 
iLaor & Dayisll20ll]) . a s oft ionizing continuu m due to intrin¬ 
sic X-ray weakness (e.g.. lLeighly et al.l2007a^ . and a soft ion¬ 
izing continuum due to small-scale absorption (e.g., Wll; 
W12). Howeyer, for the general population of PHL 1811 
analogs and WLQs, the absorption-induced soft ionizing con¬ 
tinuum appears the most likely scenario, based on systematic 
studies, albeit using small samples, of their X-ray and multi- 
wayelengthproperties (e.g., Wll; W12). 

A small-scale “shielding gas” scenario was proposed in 
Wll to explain and unify PHL 1811 analogs and WLQs, 
which was broadly motiyated by the shielding gas gener- 
ally required in the disk-wind model for BAL quasars (e.g., 
iMurrav et al.lll995l; iProga et al.ll200C)l) . IntheWl 1 scenario, 
some shielding gas interior to the BELR shields all or most 
of the BELR from the nuclear ionizing continuum, resulting 
in the obseryed weak UV emission lines. If our line of sight 
intersects the X-ray absorbing shielding gas, a PHL 1811 ana¬ 
log or an X-ray weak WLQ is obseryed; if not, an X-ray nor¬ 
mal WLQ is obseryed. Since BAL quasars were excluded 
in the selection of PHL 1811 analogs and WLQs, the incli¬ 
nation angle (with respect to disk normal) should probably 
still be relatiyely small so that the line of sight does not in¬ 
tercept the (often equatorial) disk wind which would produce 
BALs in the obseryed spectra. The reason this shielding gas 
is unusually effectiye at screening the BELR in the PHL 1811 
analogs and WLQs remains uncertain, but it should be a rare 
occurrence giyen the small numbers of PHL 1811 analogs and 
WLQs discoyered. 

The fractions of PHL 1811 analogs and WLQs among typ¬ 
ical quasars are small, < 1-2% (Wll). Howeyer, studies of 
rare and extreme objects often clearly reyeal phenomena that 
are more generally applicable, as such effe cts are more diffi - 
cult to identify in the oyerall population (cf lEddingtonll92^ . 
We shall indeed argue the case for such generality later in this 
paper (see Section 6.3 below). With the pilot studies of W11 
and W12 examining systematically the X-ray and emission¬ 
line properties of PHL 1811 analogs and WLQs, progress has 
been made toward understanding the nature of these extreme 
quasars (e.g., the Wll shielding-gas scenario). Howeyer, the 
small X-ray sample ayailable in their work, with only eight 
PHL 1811 analogs and 11 radio-quiet WLQs (which are also 
diyided into X-ray weak and X-ray normal categories for the 
case of WLQs), limited further inyestigations. A larger sam¬ 
ple is critically needed to reduce the uncertainties of stacking 
and joint spectral analyses, examine correlations between the 
degree of X-ray weakness and emission-line properties, as- 
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sess why PHL 1811 analogs are preferentially X-ray weak, 
and further explore the shielding-gas scenario. 

As an extension of the W11 and W12 work, we present here 
an X-ray and multiwavelength study of 18 PHL 1811 analogs 
and 33 WLQs, including 33 objects with new Chandra obser¬ 
vations. We describe the sample selection and Chandra data 
analysis in Sections 2 and 3, respectively. The multiwave¬ 
length properties, including the aox parameters, continuum 
SEDs, and radio properties, are presented in Section 4. In Sec¬ 
tion 5, we perform X-ray stacking and joint spectral analyses, 
estimate Eddington ratios, construct composite SDSS spectra, 
and identify spectral indicators for the X-ray weak subsample. 
In Section 6, we discuss the unification of the X-ray weak and 
X-ray normal PHL 1811 analogs and WLQs under the W11 
shielding-gas scenario, and we propose that a geometrically 
thick inner accretion disk may act as the shielding gas in the 
general population of PHL 1811 analogs and WLQs. We sum¬ 
marize in Section 7. 

We caution that, consistent with the exploratory nature of 
this work, the PHL 1811 analogs and WLQs could be a het¬ 
erogeneous population. The discussions and conclusions are 
mainly applicable to the general population of PHL 1811 
analogs and WLQs, and we acknowledge that other expla¬ 
nations are possible for a fraction of our sample. We also 
stress that PHL 1811 analogs were selected based on their 
PHL 1811-like UV emission-line properties, and they are not 
necessarily intrinsically X-ray weak like PHL 1811 itself. In 
fact, the X-ray weakness of PHL 1811 analogs is likely caused 
by absorption, based on the studies of W11 and our work 
here. Throughout this paper, we use J2000 coordinates and 
a cosmology with //n = 65.1 km s ~' Mpc“', Hm = 0.329, and 

= 0.671 (e.g.. lAde et al.ll2M5h . Pull J2000 names of the 
targets are listed in the tables while abbreviated names are 
used in the text. We quote uncertainties at a Icr confidence 
level and upper and lower limits at a 90% confidence level. 

2. SAMPLE SELECTION AND CHANDRA OBSERVATIONS 
2.1. New X-ray Sample of PHL 1811 Analogs 

Our new sample of PHL 18 11 analogs was selecte d from 
the SDSS Data Release 7 (DR7 ; lAbazaiian et al.l200^ quasar 
catalog (ISchneider et al.l 1201 011 following similar procedures 
to those in Wll, with the major difference being a re¬ 
laxation in the redshift requirement. Specifically, we re¬ 
quire redshift z > 1.7, SDSS r-band magnitude < 18.85, 
C IV A1549 REW < 10 A, C iv blueshift > 1000 km s"*, 
strong UV Pe III (UV48 A2080) and/or Fe II (2250-2650 A) 
emissionH no detection of BALs or mini-BALs (absorption 
troughs 500-2000 km s“^ wide), and radio-loudness parame¬ 
ter R < 10 (radio quiet) 0 The relaxed redshift requirement 
allows selection of a larger sample than in Wll which in¬ 
cludes more optically bright objects for efficient Chandra ob¬ 
servations. The choice of a radio-quiet sample is to avoid any 
contamination from jet-linked X-ray emission that might con¬ 
fuse the results. Some basic quasar properties (e.g., redshift, 
C IV R EW) were initial l y adop ted or derived from th e cat¬ 
alogs of iHewett & Wildl (120101) and IShen et al.l (1201 Ih . The 
C IV blueshifts were measured based on the adopted red- 

Tlie strong UV Ee II and Ee III emission relative to the other wealc UV 
lines in PHL 18 11 analogs can be exp lained in the scenario of a soft ionizing 
continuum; see lLeighlv et alH20Q7al) for details. 

The radio-loudne ss parameter is defined as R =/s GHZ //4400 A S ’ 
IKellermann et alj|1989ll . where fs gHz and /4400 A densities at 

rest-frame 5 GHz and 4400 A, respectively. 


shifts. We later performed our own meas urem ents of the 
redshifts, emission-line properties (Sec tion 12.31 below), and 
radio-loudness parameters (Section l431 below) for our Chan¬ 
dra targets. The UV Fe II and Fe III line strength was as¬ 
sessed visually among an initially selected sample satisfying 
the other criteria, and quasars with stronger Fe II and Fe III 
emission (relative to the other weak-lined quasars) were fa¬ 
vored. 

Compared to the redshift requirement in Wll,2.125 <z< 
2.385, our z > 1.7 criterion no longer requires coverage of the 
Lya and UV Fe II emission, although the visual inspection 
stage of the Fe II line strength generally selected targets with 
z < 3. An initial sample of 66 quasars was selected, including 
six of the eight radio-quiet PHL 1811 analogs in W11 (the 
other two have smaller C IV blueshifts than our criterion here). 
We ranked these objects according to the REW and blueshift 
of the C IV line, and selected 10 bright (m^ < 18.2) targets 
with the lowest C IV REWs and highest C IV blueshifts for 
Chandra observations. The selected quasars span a redshift 
range of 1.7-2.9. They were observed in Chandra Cycle 14 
with 3.7-9.5 ks exposures (Table 1) using the S3 CCD of the 
Adva nced CCD Imaging Spectrometer (ACIS: lGarmrre et alJ 

IMIH 

In addition, we obtained a 40 ks Chandra exposure of 
J15214-5202 tha t had a previous 4 ks Chandra exposure 
(lJust et al.l 120071) and was studied as a PHL 1811 analog in 
Wll. This remarkable source is one of the most luminous 
quasa rs in the SDSS catalog (M, = -30.2; e.g., lJust et al] 
120071) . and we aimed to acquire reliable basic spectroscopic 
information with this longer Chandra observation. X-ray 
properties derived from this longer observation are used in 
the relevant analyses of this study. 

2.2. New X-ray Sample of WLQs 

Our WLQ target s were mainly selected from the 
iPlotkin et al.l (1201 Obh catalog of radio-quiet WLQs which 
have REW < 5 A for all emissi on features. There wa s no C IV 
blueshift requirement for the IPlotkin et alJ (1201 Obh WLQs. 
We chose bright (/-band magnitude m, < 18.6) WLQs, ex¬ 
cluding objects that are identified as stars based on their 
proper motions, have potential absorption features (e.g., inter¬ 
vening absorption, mini-BALs, extre mely red continuum), o r 
have already been studied in Wll or[Shemmer^^Llj200^ 
There w ere 21 such WLQs selected from the IPlotkin et alJ 
(1201 Obi) catalog. In addition, we include in our WLQ 
sample two bright radio-quiet quasars found in the liter¬ 
ature that have similarly weak emission-line features, HE 
0141-3932 (iReimers et al.ll200l: REW < 15 A for C I V and 
Mg II A2799) and 2QZ J2154-3056 dLondish et all I200I 
weak [O III] A5007 and H/3 emission). These 23 WLQs span 
a redshift range of 0.5-2.5, and they were observed in Chan¬ 
dra Cycle 14 with 1.5-3.5 ks exposures using the ACIS-S3 
CCD (Table 1). 

However, of these 23 WLQ targets, J1332-(-0347 was 
later foun d to show a C IV BAL in its VLT/X-shooter 
spectrum (IPlotkin et al.j 120151) . It is also a lensed quasar 
(iMorokuma et al.l I2007^ which may have its emission-line 
REWs affected by gr avitational lensing amplification (e.g., 
IShemmer et alJ 120091) . Therefore, we will not include 

** We also observed PHL 1811 itself with a 2 Ics Chandra exposure. It 
was detected with flux and effective power-law photon index consistent with 
previous X-ray observations and did not show unexpected vaiiability (e.g., a 
return to a nominal level of X-ray emission). 
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J1332+0347 in the figures or our analyses below, since it may 
not be a bona-fide WLQ. However, we do present the basic 
properties of this object in Tables 1-3 for completeness. 

We also note that HE 0141-3932 and 2QZ J2154-3056 
were selected differently from the other WLQs. However, the 
inclusion of these two objects does not bias our results, as they 
do not have SDSS spectra and are not included in the majority 
of the statistical analyses below. Our final new X-ray sample 
of WLQs includes 22 objects. 

2.3. Redshift, UV Emission-Line, and Continuum 
Measurements 

We measured redshifts and UV emission-line properties for 
our new X-ray samples of PHL 1811 analogs and WLQs the 
same way as in Section 2.2 of Wll. Briefly, the SDSS Cat¬ 
alog Archive Server (CAS) redshifts were initially adopted. 
These redshifts may not be precise in some cases due to the 
weakness and sometimes significant blueshifts of the emis¬ 
sion lines. We examined the spectra and made small adjust¬ 
ments for some sources based on strong line features (e.g.. 
Mg II emission, narrow absorption features). We adopted 
z = 2.238 for J1521+5202 based on its best-fit H/3 line 
(Wll). The emission-line properties, including the REWs 
of the C IV A1549, Si IV A1397, A1900 complex (mainly 
C III] A1909), Ee III UV48 A2080, and Mg II A2799 emission 
features, and the blueshift and EWHM of the C IV line, were 
measured interactively. The wavelength region for each line 
was chosen based o n the upper and lower wa velength limits 
given in Table 2 of IVanden Berk el'iD (1200 Ih . and a power- 
law local continuum was fitted between the lower and upper 
10% of the wavelength region for the line measurements. 

In addition, we measured the UV Ee II REW that was not 
included in the previous studies of W11 and W12; it was mea¬ 
sured between 2250 A and 2650 A following the same ap¬ 
proach above. Eor three objects that do not have full spectral 
coverage of this wavelength range, a lower limit was derived; 
for the two cases where the covered fraction of this line fea- 
ture is > 60% based on th e SDSS quasar composite spectrum 
(IVanden Berk et aLll2001h . an estimated Ee II REW was also 
derived from this fractional coverage. The uncertainty of this 
UV Ee II REW is dominated by the continuum fitting, and a 
10% error is assumed. 

The adopted redshifts are listed in Table 1, and the UV 
emission-line properties are shown in Table 2. We also de¬ 
rived relative SDSS g-i colors, A{g-i), for our targets, which 
are defined as the co lors referenced to th e median color at 
a given redshift (e.g.,[Richa rds et al Jl2001h . The g-i colors 
were taken from ISchneider et al.l(l2010h . and the median color 
was computed within a redshift bin of Az = 0.1. The relative 
g-i colors are listed in Table 3, with positive values represent¬ 
ing relatively red spectra. We plot the individual SDSS spec¬ 
tra for the PHL 1811 analogs and WLQs i n Eigure[I] with the 
compo site spectrum of SDSS quasars from IVanden Berk et alJ 
(1200 Ih and the Hubble S p ace T elescope (HST) spectrum of 
PHL 1811 (iLeighlv et al.ll2007ah shown for comparison. 

2.4. Previous Samples Used in this Study 

In order to improve the statistics of our study, we also in¬ 
clude previous samples of radio-quiet PHL 1811 analogs and 
WLQs from W11 and W12. The basic X-ray and multiwave¬ 
length properties of these objects were adopted from Wll and 
W12. There are also X-ray sam ples of high-redshift (z » 3-6 ) 
WLQs studied previously (e.g.. lShermner et al.ll2006ll2009l) . 



Figure 1. SDSS spectra for the PHL 1811 analogs and WLQs in the Chan¬ 
dra Cycle 14 sample, following the order in Table 1 (an extended version 
of this figure showing all the spectra is available in the online journal). The 
y-axis is the flux density in arbitraiy linear units. Each spectrum has been cor¬ 
rected for Galactic extinction and smoothed using a boxcar of width 20 pixels, 
We do not show the spectrum of J1521+5202, which w as already included 
in W1 1 ■ The composite spectrum of SDSS quasars from Van den Berk et all 
<20011) and the HST spectrum of PHL 1811 ILeighlv et al 12007 al) are shown 
for comparison. 

We do not include those WLQs here due to the significant 
difference in redshift (and subsequently different rest-frame 
SDSS spectral coverage). The general comparison between 
the high-redshift WLQs and lower-redshift WLQs was pre¬ 
sented in Section 5 of W12, and their overall properties are 
consistent with each other. 

We include 7 additional radio-quiet PHL 1811 analogs 
from Wll (the other radi o-qu iet PHL 1811 analog in W11 is 
J1521+5202; see Section IZlI) and 10 radio-quiet WLQs from 
W12. There is one object, J0903+0708, in W11 that does not 
satisfy the criteria for being a PHL 1811 analog but should be 
considered a WLQ. Eor simplicity of the presentation in this 
paper, we added this object to the W12 WLQ sample. 

2.5. The Full Sample and Its Luminosity, Redshift, and 
Emission-Line Properties 

In the full sample, we have in total 18 (11 + 7) radio-quiet 
PHL 1811 analogs and 33 (22+10+1) radio-quiet WLQs. 
We show in Pigure|2]the locations of these objects in the red¬ 
shift vs. absolute /-band magnitude plane, which highlights 
the broader redshift range and brighter optical fluxes of the 
PHL 1811 analogs in our Chandra Cycle 14 sample. There is 
relatively little overlap between PHL 1811 analogs and WLQs 
in this pl ot, p ure ly d ue to the differing selection approaches 
(Sections l2.ll and lZ2]) ; physically, they are likely similar types 
of object (Wll). 

In Eigure we compare the C IV REWs and blueshifts 
of our sample objects (not all the sources in the full sample 
have C IV coverage) to those of typical radio-quiet quasars 
(iRichards et al.l l201ll) . The PHL 1811 analogs and WLQs 
have remarkably weak and strongly blueshifted C IV emis¬ 
sion, although C IV blueshift was not a selection criterion for 
the WLQs. Eor comparison, only 0.48% of the radio-quiet 
quasars in IRichards et al.l (120111) have C IV REW < 10 A, 
and 30% (1.6%) have C IV blueshift > 1000 km s“' (> 
2000 km s“'). Due to the differing selection approaches (e.g., 
Eootnote flSTl. the PHL 1811 analogs have in general some¬ 
what larger C IV REWs than the WLQs, but we do not be- 
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Figure 2. Redshift vs. absolute r-band magnitude for the samples of 
PHL 1811 analogs and WLQs. The underlying gray dots represent objects 
from the SDSS DR7 quasar catalog. The vertical dashed lines illustrate the 
nan‘o\v redshift range for the PHL 1811 analogs in W11; our ne\vly selected 
PHL 1811 analogs have a broader redshift range and generally higher fluxes 
than the Wll objects. The vertical dotted line marks the redshift (z = 1.5) 
where our sample objects start to have C IV coverage and be included in the 
C IV subsample. Solid symbols represent X-ray weak objects, while open 
symbols represent X-ray normal objects. 

lieve this indicates any fundamental difference between them; 
both groups have strikingly weak C IV emission relative to 
the general quasar population. In Figure we compare the 
Fe II REWs for the PHL 1811 analogs and WLQs. By our 
sample selection, the PHL 1811 analogs also have generally 
larger Fe II REWs than the WLQs, although there is signifi¬ 
cant overlap between the two groups. 

Since some of our WLQs are at relatively low redshifts 
{z < 1.5) with no C IV coverage in the SDSS spectra, the 
full sample might have contamination from objects having 
stronger C IV lines (e.g., iPlotkin et alJl20L5h or C IV BALs 
(e.g., similar to J1332-60347 noted above). Therefore, in some 
of the following analyses, we present results based on the sub¬ 
sample of objects at z > 1.5 with C IV coverage (the C IV sub¬ 
sample) to avoid such contamination. All the 18 PHL 1811 
analogs are in the C IV subsample, and 18 of the 33 WLQs 
(11 of the 22 in the newly observed X-ray sample of WLQs) 
are in this subsampleP^ 

3. X-RAY DATA ANALYSIS 
3.1. X-ray Counts and Photometric Properties 

We reduced and analyzed the Chandra Cycle 14 data us¬ 
ing mainly the Chandra Interactive Analysis of Observations 
(CIAO) toolsH For each source, we reprocessed the data us¬ 
ing the CHANDRA_REPRO script to apply the latest calibra¬ 
tion. The background light curve of each observation was 
inspected and background flares were removed using the DE¬ 
FLARE script with an iterative 3(7 clipping algorithm. The 
cleaned exposure times are listed in Table 1. 

For each source, we created images in the 0.5-8 keV (full), 
0.5-2 keV (soft), and 2-8 keV (hard) bands from the cleaned 
event file using the standard ASCA grade set (ASCA grades 
0, 2, 3, 4, and 6). We then ran WAVDETECT (iFreeman et all 
[20(il on the images to search for X-ray sources with a 

The C IV subsaraple also includes 11013+4927, 11604+4326, and 
12115+0001 in W12 that do not have C IV measurements in Table 2 of W12 
but actually have upper limits on the C IV REWs. 

See http : / /cxc . harvard. edu/ciao/ for details on CIAO. 


“•\/2 sequence” of wavelet scales (i.e., 1, 1.414, 2, 2.828, and 
4 pixels) and a false-positive probability threshold of 10“®. If 
a source is detected in at least one band, we adopt the WAVDE¬ 
TECT position that is closest to its SDSS position as the X-ray 
position. Seven of the 11 PHL 1811 analogs are detected, 
and 15 of the 22 WLQs are detected. The X-ray-to-optical 
positional offsets for our targets are small, ranging from 0.1" 
to 0.9" with a mean value of 0.43 ± 0.06". We also verified 
that there is no confusion with the X-ray source identification 
(e.g., no close pairs). If a source is not detected by WAVDE¬ 
TECT, the SDSS position is adopted as the X-ray position. 

We performed aperture photometry to assess the detection 
significance and extract source counts in each of the three en¬ 
ergy bands. Source counts were extracted using a 2"-radius 
circular aperture centered on the X-ray position, which corre¬ 
sponds to encircled-energy fractions (EEFs) of 0.939, 0.959, 
and 0.907 in the full, soft, and hard bands, respectively. Back¬ 
ground counts were extracted from an annular region cen¬ 
tered on the X-ray position with inner radius 10" and outer 
radius 40". For each source in each band, we computed a 
binomial no-source proba bility, Pr, to asses s the significance 
of the source si gnal (e.g.. lBroos et al.ir2007t IXue et al.ll201 lb 
iLuo et aO2013l) . which is defined as 


Pb{X>S) = J2 


m 


x=s 


x\{N-xy. 


p^il-p) 


,N-X 


( 1 ) 


In this equation, S is the total number of counts in the source- 
extraction region; N = S-\- B, where B is the total number 
of counts in the background extraction region; p = 1/(1 -I- 
BACKSCAL), where BACKSCAL is the ratio of the area of the 
background region to that of the source region. Pb represents 
the probability of observing the source counts by chance un¬ 
der the assumption that there is no real source at the relevant 
location. If the Pb value in a band is smaller than 0.01 (cor¬ 
responding to a > 2.6(7 detection), we considered the source 
detected in this band and calculated the Icr errors on the net 
counts, which were derived from the Icr errors on the ex¬ 
tracted source and background counts (iGehrels [1^8^ follow¬ 
ing the numerical method in Section 1.7.3 of iLvonsI (119911) . 
We note that this Pb threshold is appropriate in our case where 
the position of interest is precisely specified in advance. If the 
Pb value is larger than 0.01, we considered the source unde¬ 
tected and derived an up per limit on the sou rce counts using 
the Bayesian approach of iKraft et al.l (119911) for a 90% confi¬ 
dence level. In Table 1, we list the source counts and upper 
limits in the soft and hard bands. 

We derived a 0.5-8 keV effective power-law photon in¬ 
dex (Feff) for each source from the band ratio between the 
hard and soft band counts, calibrated using the Portable, In¬ 
teractive, Multi-Mission Simulator (PIMMSiFI with the as¬ 
sumption ofa_gower4aw_s2ectrum^odifled with Galactic ab¬ 
sorption (iDickev & Lockm^ll990l) . Feff represents the real 
power-law photon index if the underlying spectrum is truly 
a power law with Galactic absorption. The uncertainties of 
(or limits on) the band ratios (and subs equently F^ff) we re de¬ 
rived using the Bayesian code BEHR (iPark et alJl2006h . For 
sources undetected in both the soft and hard bands, F^f can¬ 
not be constrained and a value of 1.9 wa s ad opted for flux 
or flux density computations later ('Section lTTI b elow), which 
is typical for luminous radio-quiet quasars (e.g.- lReeves et alJ 
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Figures, (a) C IV REW vs. C IV blueshift for the samples of PHL 1811 analogs an d WLQs. Only object s having both C IV REW and blueshift measurements 
are shown here. The underlying gray dots are the 13 582 typical r adio- qui et qu asars in IRichards et alj j2Ql ID . The PHL 1811 analogs and WLQs have remarkably 
weak (as expected from the selection criteria detailed in Sections [2.1l and [2^ and usually strongly blueshifted C IV emission, (b) C IV REW vs. Fe II REW for 
the samples of PHL 1811 analogs and WLQs. An'ows represent upper limits (for C IV REWs) or lower limits (for Fe II REWs). The PHL 1811 analogs have in 
general larger Fe II REWs than the WLQs by selection. In both panels, solid symbols represent X-ray weak objects, while open symbols represent X-ray normal 
objects; the green star represents PHL 1811. 


119971 Uust et al.ll2007l : IScott et al.ll201 Ih . The band ratios and 
Teff values are listed in Table 1. 

3.2. Spectral Analysis of J1521+5202 

J1521+5202 is an exceptionally luminous quasar (e.g., Fig¬ 
ure la. However, it is remar kably X-ray weak (by a factor 
of Ri 35), as first reported by Llust et all (120071) using a 4 ks 
Chandra observation. As a PHL 1811 analog, J1521-I-5202 
was studied in Wll, which also presented a near-infrared 
(NIR) spectrum showing very weak [O III] emission and a 
higher optical Fe Il/H^ ratio than those of typical quasars. 
The source was only weakly detected in the 4 ks Chandra 
observation with 3 full-band counts, preventing further con¬ 
straints on the band ratio and effective photon index. With 
our much longer Chandra observation (37.1 ks cleaned expo¬ 
sure) of J1521-1-5202, we derived a band ratio of 1. Vg 2 and an 
effective power-law photon index of 0.6 ± 0.2 (Table 1). This 
flat effective photon index, compared to the ty pical value of 
F R! 1.9 for luminous radio-quiet quas ars re.g.. lReeves et alJ 
[TW^LTiist et al .11200^ l^t et al 1201 Ih . suggests that the sig¬ 
nificant X-ray weakness of this source is likely due to strong 
absorption. 

We then performed basic spectral analysis for J1521-I-5202. 
The 0.3-8 keV (rest-frame 1.0-25.5 keV) spectrum was ex¬ 
tracted using the CIAO SPECEXTRACT script within a circular 
aperture of 4" radius centered on the X-ray position. There 
are r 90 source counts extracted including < 1 background 
count. The background spectrum was extracted from an annu¬ 
lar region with inner radius 6" and outer radius 15", which is 
free of X-ray source contamination. Given the limited source 
counts, we fit ted the spectrum using XSPEC (version 12.8.1; 
lArnaudiri996l) with a simple model consisting of a power law 
modified by both intrinsic absorption (at z = 2.24) and Galac¬ 
tic absorption. The C statistic (c stat) was used due to the lim¬ 
ited photon counts (ICashlll97^ . The best-fit photon index is 
F = 1.5Vg 35 , and the intrinsic absorption column density is 
Ah = (1.31];o4f) X 10^^ cm“^ (C = 88 for 75 degrees of free¬ 
dom). The spectrum and the best-fit model are displayed in 
Figure |4] There is a possible emission-line feature at rest- 
frame 7.5 keV, which wou ld likely correspond to a blueshifted 
Fe line if confirmed (e.g.. [Reeves et al.ll2004^ . More data are 



Observed-Frame Energy (keV) 


Figure 4. Chandra spectrum of J1521+5202 (blue) overlaid with the best-fit 
model (red). The 0.3-8 IceV spectrum was fitted with an absorbed power-law 
model modified with Galactic absorption using cstat in XSPEC. For display 
purposes, the data are grouped with a minimum of six counts per bin. The 
inset shows the best-fit F and N-n values (magenta plus), and their la (blue), 
2a (green), and 3a (red) contours. The hottom panel shows the ratio of the 
data to the hest-fit model. The hest-fit model suggests strong intrinsic X-ray 
absorption, with Nh = (1.3±0.5)x 10“^ cm“^. 

required to constrain this feature. 

The spectral analysis also suggests significant X-ray ab¬ 
sorption in J1521+5202, in agreement with the simple flat 
Feff argument above. We note that the absorption derived 
from our simple spectral model cannot entirely account for 
the observed X-ray weakness; after correction for this absorp¬ 
tion J1521+5202 is still a factor of R 7 times X-ray weak (a 
R 2.5(7 deviation from the aox-Tjjgg ^ relation). This result 
suggests that our simple model probably did not recover the 
real (larger) absorption column density associated with a more 
complex absorber, as is commonly seen in highly obscured 
AGNs 13 Alternatively, j 1521+5202 may also be intrinsically 

Therefore, the best-fit F value is liJcely not representative of the intrinsic 
power-law spectral shape, and we cannot obtain a reliable estim ate of the 
Eddington ratio based on this X-ray photon index (see Section [531 for further 
discussion). 
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X-ray weak by a factor of a few besides the strong X-ray ab¬ 
sorption. 

4. MULTIWAVELENGTH PROPERTIES 
4.1. Distribution of the X-ray-to-Optical Power-Law Slopes 

We measured the aox parameters for our Chandra Cy¬ 
cle 14 sample objects. Th e rest-frame 2500 A flux densi¬ 
ties were adopted from the IShen et aTI (1201 ll) SDSS quasar 
catalog; for three objects that were not included in this cat¬ 
alog, the flux densities were interpolated/extrapolated from 
their optical photometric data. The rest-frame 2 keV flux den¬ 
sities were derived from the soft-band net counts (Table 1) 
using PIMMS assuming a power-law spectrum modified with 
Galactic absorption; the observed soft band (0.5-2 keV) cov¬ 
ers rest-frame 2 keV for our sources. If a source is not de¬ 
tected in the soft band, an upper limit on /2 keV was calcu¬ 
lated. There are two sources, J0147H-0003 and J1133H-1142, 
that are not detected in the soft band but are detected in the 
hard band. We still present upper limits on /2 keV (and aox) 
for these two sources, although the /2 keV flux density could 
be estimated via extrapolation of the hard-band flux assum¬ 
ing a power-law photon index. This approach does not af¬ 
fect the analyses in this study. The aox values, along with 
other relevant parameters, are listed in Table 3. The level 
of X-ray weakness is usually measured by the Aaox pa¬ 
rameter, which is defined as the difference between the ob¬ 
served aox and the one expected from the aox-T 250 () A 
tion (Aaox = aox.obs-aox.exp)- From Aaox, we can derive 
the factor of X-ray weakness (/weak = 403“'^“™) with respect 
to the expected X-ray flux. These parameters are also listed in 
Table 3. 

The aox vs. Tpano A distribution for our full sample is dis¬ 
played in Figure |5] also shown for comparison are PHL 1811 
and some typical AGN sam ples (ISteffen et al.l200(itl.lust et ^ 
120071; iGibson et al.ll2008h . We adopt Aaox = -0.2 (/weak = 
3.3) to be the dividing threshold between X-ray weak and X- 
ray normal quasars, which corresponds to a Ri 1.3 ct (~ 90% 
single-sided con fidence level) offset given the aox rms scat¬ 
ter in Table 5 of ISteffen et alJ (|2006|) . A remarkable fraction 
of our sample sources are X-ray weak, similar to PHL 1811. 
Specifically, of the 18 PHL 1811 analogs, 17 (94!^3%) areX- 
ray weak, and 16 (48);[2%) of the 33 WLQs are X-ray weak. 
In the C IV subsample, 8 (44:^ij%) of the 18 WLQs are X-ray 
weak. For comparison, t he fraction of X-ray weak objects in 
the improved version of iGibson et al.l (l2008ll Sample B (see 
Footnote 16 of Wll) is only 7.6% (10/132). 

The significantly different fractions of X-ray weak objects 
among typical quasars and our samples of PHL 1811 analogs 
and WLQs are also evident in the Aaox distributions (Fig¬ 
ure 01. As in Wll and W12, we used the Peto-Prentice 
test implemented i n the Astronomy Survival Analysis pack¬ 
age (A SURV; e.g.. iFeigelson & NelsonllT9^ iLayalley et al.l 
Il992h to assess whether two samples follow the same distri¬ 
bution. For PHL 1811 analogs, the probability of their Aaox 
yalues being drawn from the same parent population as those 
of typical quasars (improyed yersion of IGibson et al.l l2008l 
Sample B) is Ri 10“^^ (10.4cr). For WLQs, this probability 
is 3 X 10-5 (4.1ct). 

Since most of our X-ray weak PHL 1811 analogs and 
WLQs are not X-ray detected and only upper limits on Aaox 
are ayai lable , we also performed X-ray stacking analyses (see 
Section lSTI below) to obtain the ayerage leyel of X-ray weak¬ 
ness. The stacking results are listed in Table 4 as subsamples 



29.0 29.5 30.0 30.5 31.0 31.5 32.0 32.5 

log L(2500A) (ergs s"' Hz"') 


Figure 5. X-ray-to-optical power-law slope (aox) vs. 2500 A monochro¬ 
matic luminosity for the samples of PHL 1811 analogs and WLQs. For 
compar ison, the 132 radio- quiet, non-BAL quasars in the improved ver¬ 
sion of IGibson et alj 12008ft Sample B (see Eootnote 16 of Wll) are rep¬ 
resented by the black dots, and PHL 1811 is shown as the green star. The 
s mall grey dots and do wnw ard arrow s (uppe r limits) are from the samples 
of ISteffen et all 42006ft and lJust et al.l 12007ft . and the solid line represents 
th mJusTet alj 120071) anY-L 2 ;oo a fElation- Th® dashed line (Aaox = -0.2) 
marks the division between X-ray weak and X-ray normal quasars adopted in 
this study. 

“P1811A SB undet” and “WLQ undet” for the undetected 
sources (excluding two XMM-Newton undetected sources), 
and the stacked Aaox yalues are plotted as gray bars in Fig¬ 
ure 0 The corresponding ayerage X-ray weakness factors for 
the undetected PHL 1811 analogs and WLQs are 513 and 108, 
respectiyely. As noted aboye, there are two PHL 1811 analogs 
in the stacked sample, J0147-I-0003 and J1133-1-1142, that are 
not detected in the soft band but are detected in the hard band. 
Excluding these two sources from the stacking still yields a 
large negatiye stacked Aaox yalue (-0.87), as shown in the 
results for subsample “P1811A undet” in Table 4. 

Including detected objects, the stacked Aaox yalues are 
-0.61 and -0.47 for the X-ray weak PHL 1811 analogs and 
WLQs, corresponding to ayerage X-ray weakness factors of 
39 and 17, respectiyely. The X-ray weak WLQs in the C IV 
subsample have a similar stacked Aaox value (-0.62; average 
X-ray weakness factor of Ri 41) to the X-ray weak PHL 1811 
analogs. In an absorption scenario, an X-ray weakness factor 
of 40 corresponds to Ah ~ 9 x 10 ^^ ^^ 1-2 fop ^ = 2 and F = 1.9 
in the MYTORUS model (iMurphv & Yaaoobll2009h . Using 
the Kaplan-Meier estimator in ASURV, we also derived the 
mean Aaox values for the X-ray weak PHL 1811 analogs and 
WLQs, which are -0.49 ±0.03 and-0.38±0.02, respectively. 
Note that the Kaplan-Meier estimator assumes that the cen- 
sored data have an identica l distribution as the measured data 
(IFeigelson & Nels^ll985l) . which may not be applicable to 
our Aaox values here; also, the mean and stacked values dif¬ 
fer in the sense that the mean Aaox estimates the geometric 
mean of the X-ray weakness factors while the stacked Aaox 
estimates the arithmetic mean. 


4.2. IR-to-X-ray Spectral Energy Distributions 

We constructed infrared (IR) to X-ray continuum SEDs 
for our full sample of PHL 1811 analogs and WLQs, us¬ 
ing photometric data c ollected from the W ide-field Infrared 
Survey Explorer ('W/5£ ’;IWright et £0120 lOl). Two Micron All 
Sky Survey (2MASS; ISkrutskie et al, 20061). SP SS, and/or 
Galaxy Evolution Explorer (GALEX; iMartin et al.ll2005h cat- 
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Figure 6. Distribution of the Aoqx values for the (a) PHL 1811 analogs and (b) WLQs in our sample. The solid and hatched shaded histograms (thick and 
thin leftward arrows) represent X-ray detected (undetected) sources in the Chandra Cycle 14 sample and the previous W1 1 or W12 sample, respectively. In both 
panels, the gray bar shows the average Aoqx value for the undetected sources estimated via stacking analyses (see Sectio nIS.llbelow): there is one XMM-Newton 
undetected source that was not included in the stacking in each case. For comparison, the Aoqx distribution for the 132 IGibson et all <20081) Sample B quasars 
is presented as the unshaded histogram, and the Aaox value of PHL 1811 is shown as the green star. A significant fraction of the PHL 1811 analogs and WLQs 
are X-ray weak compared to typical quasars. 


alogs. The optical and UV data have been corrected for 
Galactic extinction follo wing the deredden ing approach of 
ICardelli et al.l (Il989l) and lO’Donnelll (11994) . The combined 
SEDs for the PHL 1811 analogs and WLQs are displayed 
in Ligures |7] and |8] respectively. The mean SEP of op ti- 
cally luminous SDSS quasars from [Richards et all (120061) is 
shown for comparison. Because of our sample selection cri¬ 
teria, the PHL 1811 analogs are about a fact or of fti 3 times 
more optically luminous on average than the iRichards et al.l 
( 120061) optically luminous SDSS quasars, while the WLQs 
have comparable luminosities to the SDSS quasars. Regard¬ 
less of whether they are X-ray weak or X-ray normal, our 
PHL 1811 analogs and WLQs in general have IR-to-UV SEDs 
similar to those of typical quasars. Such SEDs differ signifi¬ 
cantly from those of BL Lac objects, indicating that our sam¬ 
ple is not contaminated by any radio-quiet BL Lac candidates 
(e.g., see discussion on J1109-)-3736 in W12). Studies of a 
sample of high-r edshift (z = 2.7-5 .9) WLQs also found typi¬ 
cal quasar SEDs (iLane et al .11201 ll) . 

We estimated bolometric l uminosities (Lbui) for our sam¬ 
ple objects by integrating the IRichards et al.l (l2006h compos¬ 
ite SED (between lO'^'^-lO'® Hz) scaled to their rest-frame 
3000 A luminosities^ The resulting bolometric luminosi¬ 
ties for the full sample are listed in Table 5, ranging from 
Ri 2 X 10^^ erg s“' to Ri 2 x 10"^* erg s“' with a median value 
of Ri 8 X 10"^® erg s“'. We compared these bolometric lumi¬ 
nosities to those estimated via integrating the SEDs directly. 
The ratios of the two luminosities have a median value of 
Ri 1, and an rms scatter of Ri 3, suggesting that the typical 
uncertainty associated with the estimated LboI values is ap¬ 
proximately a factor of three. There may also be additional 
systematic uncertainty if the EUV properties of these obj ects 
differ from those of typica l quasars (e.g., see Sections [6. H and 
16.21 iKrawczvk et al.ll2013h . 

4.3. Radio Properties 


The 0.5-1 0 keV luminosi t y only constitutes 4.5% of the bolometric lu¬ 
minosity in the (Richards et all <2006h composite SED. Therefore, assuming 
a nominal level of X-ray emission for our sample objects when estimating 
bolometric luminosities does not affect the results materially. 
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Figure?. Combined SED for the PHL 1811 analogs. The IR-to-UV SED 
data were collected from the WISE (red), 2MASS (magenta), SDSS (blue), 
and GALEX (orange) catalogs. These data have been corrected for Galac- 
tic extinction f o llowin g the dereddening approach of ICardelli et al.l J1989I) 
and tO’Donnelll 419941) . The green data points and arrows show the 2 keV 
luminosities and upper limits. The SED for each obj ect was scaled to the 
composite quasar SED of optically luminous quasars 4Richards et al.|[200^ 
at rest-frame 3000 A, and then combined. Some of the GALEX data points 
affected by the Lyman break are removed in the combined SED for display 
purposes. The combined IR-to-UV SED of PHL 1811 analogs is similai* to 
the composite quasar SED. We do not separate the X-ray normal PHL 1811 
analog (J1537+2716) from the other X-ray weak objects as it also has a typi¬ 
cal quasar SED. 


Combine(j SED of 18 : 



2 keV ♦ J 


- I-Richards et al. (2006) quasar SED | '' - 

..I.I I.ll I I' 

3 14 15 16 17 1 

log i^rest (Hz) 


Only radio-quiet {R < 10) quasars are included in our sam¬ 
ple. The radio-lo udness parameter was c omputed as R = 
/s GHz //4400 A le.g.. lKellermann et al.lll989h . where the rest- 
frame 5 GHz and 4400 A flux densities were converted from 
the observed 1.4 GHz and rest-frame 2500 A flux densities 
(or their upper limits) ass uming a radio power-law slope of 
Qr = - 0.8 oc e.g., iLalcke et alJll996l: iBarvainis et alJ 
20051) and an optical p ower-law slope of = “0.5 (e.g., 
Vanden Berk et al.l 120011) . We obtained radio flux informa¬ 
tion at 1.4 GHz via cross-matching to the Laint Images of 
the Radio Sky at Twen t y-Centimeters (LI RST) survey cata¬ 
log (lBeckeFet^|T995|; I^MtEiESJlIlO^). Lor sources not 
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Figures. Same as Figure [7] bu t for th e (a) X-ray weak and (b) X-ray normal WLQs. The black stars represent the composite SED of 18 high-redshift 
(z = 2.7-5.9) WLQs in lLane et alj <201111 . The combined SEDs of the X-ray weak and X-ray normal WLQs are similar from the IR to UV, and they agree well 
with the composite quasar SED. 


detected by the FIRST survey, the upper limits on the radio 
fluxes were set to 0.25-t3tTi-ms mJy, where arms is the rms noise 
of the FIRST survey at the so urce position and 0 .25 mJy is to 
account for the CLEAN bias (IWhite et al.lll997h . The R val¬ 
ues for the Chandra Cycle 14 targets are listed in Table 3. 

Given the radio-quiet selection criterion, all our sample 
objects have optical-to-radio power-law slopes ctro > -0.21 
(e.g.. Equation 4 of W12) and are in the radio-quiet region of 
the aox vs. ctro plot (Eigure 7 of W12). However, the two 
objects with the highest R values in our sample (J0844H-1245 
with R = 8.0 and J1156H-1848 with R = 7.4; see Table 3) 
are considered radio-intermediate in IShen et all (1201 ll) with 
f^shen = 11.2 and Rshen = 10.6, respectively. The radio- 
loudne ss parameter was defined differently in IShen et aTI 
( 1201 Ih with Rshen = /s ghz//25oo A’ which leads to the above 
small discrepancy. Both quasars are X-ray normal, and we 
cannot exclude entirely the effect of jet contamination in these 
two objects. In fact, J0844H-1245 is an outlier in the diagnostic 
plot for X-ray weak quasars (Eigure fTblbelow'). which might 
be related to its relatively high radio loudness. Since both 
objects are still radio-quiet under our definition, and we have 
also verified that the inclusion of the two does not affect our 
statistical analyses significantly, we do not remove them from 
our study. 

5. SAMPLE ANALYSES 

5.1. X-ray Spectral Stacking Analyses of the X-ray Weak 
Objects 

The X-ray spectral analysis of J1521H-5202 (Section l3.2t 
suggested significant X-ray absorption. For the other X-ray 
weak objects in our sample, there are few X-ray photons de¬ 
tected for spectral analyses. The derived effective photon in¬ 
dices (Table 1) also have large uncertainties and cannot con¬ 
strain whether a source has a flat X-ray spectrum that is in¬ 
dicative of absorption. In this case. X-ray stacking analyses 
are often used to assess the average spectral properties of the 
sample. With this approach, Wll and W12 found relatively 
flat/hard stacked effective photon indices for their samples, al¬ 
though with large error bars, suggestive of X-ray absorption. 
With the addition of our new X-ray samples, we can improve 
the statistics of the stacking results and obtain tighter spectral 
constraints. 

We stacked the X-ray photometry by combining the ex¬ 


tracted source and background counts of indiv idua l sources 
and following the same procedure as in Section [TTI to derive 
photometric properties of the stacked source. Since our ob¬ 
jects are generally at high redshifts (a mean value of z = 1.85 
for the full X-ray weak sample), we are probing the highly 
penetrating X-ray emission in the Ri 2-23 keV rest-frame 
band that is more reliable for assessing the presence of heavy 
absorption than lower-energy observations. Stacking analy¬ 
ses were performed for several subsamples of X-ray weak 
objects; the results are given in Table 4, which include the 
number of sources used in the stacking, mean redshift, total 
stacked exposure time, stacked soft- and hard-band counts, 
stacked effective photon index, and stacked aox and Aaox- 
For the full X-ray weak sample, the stacked effective photon 
index is relatively hard (Feff = 1 .37^g 23 ) compared to the typi- 
cal value of F ~ 1 . 9 for luminous r adio-quiet quasar s (e.g., 
iReeves et alJll997l Llust et akl 120071 : IScott et aklboi Ih . The 
hard average spectral shape is more evident in the subsam¬ 
ple of X-ray weak PHL 1811 analogs, with Feff = 

The stacking results suggest that, unlike PHL 1811 itself that 
appears intrinsically X-ray weak (at least in the 0.5-8 keV 
band), these X-ray weak PHL 1811 analogs and WLQs on av¬ 
erage are X-ray absorbed, consistent with the con clusi on from 
the spectral analysis of J1521H-5202 in Section [LT] How¬ 
ever, given the uncertainties associated with the stacking pro¬ 
cedure (e.g., individual sources contribute differently to the 
stacked signal due to the different fluxes, exposure times, and 
rest-frame bands probed), we cannot exclude the possibility 
that a fraction of the sources are intrinsically X-ray weak like 
PHL 1811 (e.g., those sources with no or poor constraints on 
Left). 

5.2. Joint Spectral Fitting of the X-ray Normal Objects 

The X-ray stacking analyses of the X-ray weak objects 
likely did not provide constraints on their intrinsic X-ray spec¬ 
tral shapes (due to the likely modification of their spectra by 
X-ray absorption), and these objects have too few X-ray pho¬ 
tons for spectral fitting (the stacked number of counts for all 
the X-ray weak objects is only 87). Therefore, to obtain a rel¬ 
atively tight constraint on the intrinsic X-ray power-law pho¬ 
ton indices of PHL 1811 analogs and WLQs, we performed 
joint spectral fitting of the one X-ray normal PHL 1811 analog 
and 17 X-ray normal WLQs. The individual source and back- 
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Figure 9. Contour plot for the photon index and intrinsic absorption col¬ 
umn density from the joint spectral analysis for the X-ray normal PHL 1811 
analogs and WLQs. The blue, green, and red curves are the Icr, 2(j, 3cr con¬ 
tours, respectively. The magenta plus sign indicates the best-fit F and A^h 
values. 

ground spectr a we re extracted following the same procedure 
as in Section 13.21 and all eighteen spectra in the rest-frame 
> 2 keV band were htted jointly with an absorbed power-law 
model using cstat in XSPEC, allowing each source to have its 
own redshift and Galactic column density. There are Ri 610 
total spectral counts for the htting including only r; 2 back¬ 
ground counts, and these spectral counts are not dominated 
by the counts from one or a few objects. 

The best-fit photon index is T = 2.18 ±0.09 (C = 449 for 
487 degrees of freedom), and no intrinsic absorption is re¬ 
quired with an upper limit of A^h < 2.9 x 10^' cm“^. The 
1-3(7 contours for the best-fit T and A^h values are shown 
in Figure |9] The same analysis was also performed for the 
C IV subsample, which includes one X-ray normal PHL 1811 
analog and nine X-ray normal WLQs. The resulting photon 
index is slightly higher, T = 2.26 ±0.11 (C = 283 for 330 
degrees of freedom), and no intrinsic absorption is required 
(A^h < 1.4 X 10^' cm“^). Similar joint fitting was performed 
for sev en high-redshift radio-quiet WLQs in IShemmer et al.l 
(l2009ll . and their resulting L constraint is consistent with our 
r values within the errors, although our uncertainties are fac¬ 
tors of Ri 5 times smaller owing to the substantially larger 
number of spectral counts available. 

5.3. Eddington-Ratio Estimates 

The extreme emission-line properties of PHL 1811 
analogs and WLQs could perhaps be related to rapid or 
even super-Eddington accretion (e.g., iLeighlv et al.ll2007allbt 
IShemmer et al.l 1201 Ot Wll), and thus we estimated their 
Eddington ratios (LBoi/TEdcj) based on their estimated bolo- 
metric luminosities and supermassive black hole (SMBH) 
masses. Most of our PHL 1811 a nalogs and WLQs have 
their SMBH masses estimated in IShen et aT1 (1201 ll) us¬ 
ing the single-epoch virial mass approach, with 25 ob¬ 
jects having Mg ll-based estimates and 21 objects hav¬ 
ing C iv-based estimates. The C IV virial mass approach 
is likely not applicable for our objects as the prominent 
C IV blueshifts i ndicate a strong wind component for the 
C IV BELR (e.g..lBaskin & LaoifcoOSHRichards et akllMlt 

iTrakhtenbrot & Netzeii 120121: IShenI I 2 OI 3 I) . and C iv-based 

masses may be an order of magnitude hig her than the real val- 
ues for quasars with weak C IV lines ('e.g.. lKratzer &. RichardsI 


IMl- 

For the 25 objects with Mg ll-based SMBH masses, the 
median LBoi/TEdd value is 0.27 with an interquartile range 
of 0.17 to 0.38 (the mean Luoi/TEdd value is 0.39). Even 
the Mg ll-based SMBH masses are perhaps not reliable for 
our objects that have weak Mg II line emission, although 
not as abnormally weak as the C IV line emission (e.g., 
iPlotkin et al.l I 20 1 5h 1^ A recent VLT/X-shooter study of 
WLQs dPlotkin et al.ll2015[) also suggests that, at least some¬ 
times, the Mg II BELR is complex and may not be virialized 
in these exceptional objects. 

Four of the objects with Mg ll-based masses (J0945H-1009, 
J1411H-1402, J1417H-0733, and J1447-0203) recently have 
had their masses estimated using the H/3 line profiles 
(IPlotkin et al.l2015h . and the masses are on average smaller by 
a factor of r: 3. The H/3 lines of WLQs are generally weak at 
a similar level to their Mg II lines (IPlotkin et al.l2015h , but the 
H/3 mass estimator is considered the most reliable one among 
all the single-epoch virial mass approaches (e.g.. lShenll2013[ 
and references therein). The estimated Eddington ratios for 
these four objects plus J1521H-5202 which also has an H/3- 
based mass (Wll) range from 0.9 to 2.8 with a median value 
of 1.4 (a mean value of 1. 7)0 The available SMBH mass 
and Eddington-ratio estimates for our full sample are listed 
in Table 5. The uncertainties on th ese SMBH masse s include 
measurement errors (Ri 0.15 dex; IShen et aT]l201 ll) . system- 
atic errors of the virial-mass approach (r; 0.4-0.5 dex; e.g., 
IShenllMT^ . and additional systematic errors (unknown but 
likely large) related to the likely unusual BELRs of our ex¬ 
treme quasars. Given the substantial uncertainties associated 
with the SMBH masses and likel y also the estimated bolomet- 
ric luminosities (see Section l±2l l. the estimated Eddington ra¬ 
tios do not provide strong constraints on the accretion status 
of the PHL 1811 analogs and WLQs. 

We also make basic estimates of Lnoi/TEdd from the hard 
X-ray power-law phot on in dices obtained from our joint spec¬ 
tral analyses (Section I5.21 i. which have been demonstrated 
to be a useful direc t probe of the Eddi n gton ratio (e.g., 
Shemmer et al.l l2008t [Risaliti et al.l l2009t Brightman et al.l 
2 OI 3 I) . Quasars with high Eddington ratios generally have 
very soft X-ray photon indices. For example, PHL 181 1 
has r = 2.3 ±0.1 and Luoi/TEdd ~ 1.6 (ILeighlv et al.ll20 07bh. 
Based on the empirical L-LboI /Tpdd relations (IShemmer et al.l 
I 2 OO 8 I; iRisaliti et al.l 1200^ iBrightman et al.ll2013l) . the best- 
ht hard X-ray photon index for our X-ray normal objects, 
r R 2.2, corresponds to Luoi/TEdd ~ 1- Therefore, the joint 
spectral fitting results suggest high Eddington ratios for our 
X-ray normal objects, or our PHL 1811 analogs and WLQs 
in general considering that the X-ray weak and X-ray nor¬ 
mal objects could be unihed (Wll; see also Section 6 be¬ 
low). There is substantial object-to-object intrinsic scatter in 
the T-LBoi/TEdd relations. Thus, by using the average T from 
the joint spectral fitting of 10-18 objects, our constraint on the 
typical LboI /^Edd should be more reliable than that from the T 
values of a few individual objects. 

In summary, the Eddington-ratio estimates derived from the 
virial SMBH masses are highly uncertain and perhaps sys- 

Being a high-ionization line, the C IV line of PHL 1811 analogs 
and WLQs is expected to be more significantly affected by a soft ionizing 
SEP compared to lo w-ionization lines such as Mg II (e.g., Section 4.1.4 of 
ILeighlv et al.ll2007aft . 

T hese Eddington ratios differ from those in Wll and IPlotkin et al.l 
<201511 as the bolometric luminosities were estimated via different ap¬ 
proaches. 
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tematically in error. We caution that individual measurements 
should not be over-interpreted and, at best, one can look at 
the group properties for some average indication. The group 
properties from the H/3-based virial masses suggest the Ed¬ 
dington ratio may be high. Meanwhile, the large hard X-ray 
photon index from the joint spectral fitting suggests a high 
Eddington ratio for our objects as a group. We consider it 
plausible that our PHL 1811 analogs and WLQs are accreting 
at high or even super-Eddington rates. 

5 . 4 . Composite SDSS Spectra 

We constructed composite SDSS spectra for our PHL 1811 
analogs and WLQs, to compare their average UV spectral 
properties to those of PHL 1811 and typical quasars and also 
to make comparisons between the X-ray weak and X-ray nor¬ 
mal populations within our sample. Median composi te spec¬ 
tra were created following IVanden Berk et all (1200 Ih . Basi¬ 
cally, each individual spectrum was corrected for Galactic ex¬ 
tinction, smoothed using a boxcar of width 10 pixels, and nor¬ 
malized at rest-frame 2240 A. The continuum at rest-frame 
2240 A is not significantly co ntaminated by Ee II or other 
lines (e.g.- IRafiee & Hallll2011h . The composite spectra were 
derived by computing the median flux density at each wave¬ 
length where at least five objects have spectral coverage. 

The composite SDSS spectrum for the 17 X-ray 
weak PHL 1811 analog s is shown in Eigure [TOk : the 
IVanden Berk et al] (1200 lb composite spectr um of SPSS 
quasar s and the HST spectrum of PHL 1811 (iLeighlv et alJ 
l2007ab are included for comparison. The weak emission lines, 
including C IV, C III] , and Mg II, are evident in PHL 1811 and 
its analogs. Meanwhile, the PHL 1811 analogs have similar 
Le III UV48 A2080 and UV Le II (2250-2650 A) emission 
to the SDSS composite spectrum, while PHL 1811 has even 
stronger UV Le III and Le II emission. On average, PHL 1811 
and its analogs have similar continua that are redder than the 
SDSS composite spectrum between r; 1500-2700 a 13 Us¬ 
ing the SDSS filter informatior^] and the median redshift of 
z = 2.2, we measured a relative g-i color (in the observed 
frame) of A(g- /) = 0.24 for the PHL 1811 analog composite 
spectrum. 

To check the reliability of our approach for creating com¬ 
posite spectra, we made a composite spectrum for the 7251 
SDSS quasars satisfying the redshift, magnitude, radio¬ 
loudness, and non-BAL criteria of PHL 1811 analogs (Sec¬ 
tion 12. lb . and it agrees very well with the IVanden Berk et al.l 
(120011) SDSS quasar composite spectrum in the 1500-3500 A 
range. In addition, to check whether the redder color [A(g- 
i) = 0.24] of the 17 X-ray weak PHL 1811 analogs is an out¬ 
come of small sample bias, we randomly picked 17 SDSS 
quasars from the above sample and measured the A(g - i) 
value of their composite spectrum. This practice was repeated 
10 000 times and the resulting A{g-i) distribution has an rms 
scatter of 0.037. Therefore, the redder color of the PHL 1811 
analogs is a highly significant (6Aa) result. 

The composite SDSS spectra for the 14 X-ray weak and 
17 X-ray normal WLQs are shown in Eigure [TOb. with com¬ 
parison to the SDSS quasar composite spectrum. Similar to 
the PHL 1811 analogs, the X-ray weak WLQs have weak 


C IV, C III], and Mg II emission lines but normal UV Le III 
and Le II emission, compared to the SDSS composite spec¬ 
trum. The X-ray weak WLQs have generally redder con¬ 
tinua than the SDSS composite spectrum, with A(g-/) = 0.11. 
The X-ray normal WLQ composite spectrum is bluer than the 
SDSS composite spectrum, with A(g — i) = -0.1 10 It also 
has weaker UV Le II emission than the X-ray weak WLQ or 
SDSS composite spectra (more evident in Ligures fT^ and fOl 
below). 

Composite spectra were also made for the C IV subsample. 
In Eigure [TT] we compare the composite spectra for the X-ray 
weak PHL 1811 analogs. X-ray weak WLQs, and X-ray nor¬ 
mal WLQs in the C IV subsample. These spectra are similar to 
the corresponding composite spectra for the full sample. The 
composite spectra for the X-ray weak PHL 1811 analogs and 
WLQs are similar to each other0 The X-ray weak WLQs 
have redder continua and stronger UV Le II (2250-2650 A) 
emission than the X-ray normal WLQs. 

We de-reddened the composite spectrum for the X-ray weak 
WLQs with various Ay values to test whether the difference 
between the X-ray weak and X-ray normal composite spec¬ 
tra can be explained by dust extinction. A Small Magellanic 
Cloud (SMC) extinction law (iGordon et al.ll2003l : Ry = 2.74) 
was adopted w hich is usually used to describe intrinsic qu asar 
reddening (e.g.. [Hopkins et al.l2004tlGlikman et al.l2013 . As 
shown in Eigure [12] mild dust extinction of Ay = 0.05 mag 
could describe well the redder composite spectrum for the 
X-ray weak WLQs, with the remaining residuals mainly 
in the UV Le II and Mg II emission. This small amount 
of dust reddening corresponds to a gas column density of 
Nh < 10^° cm“^ assuming a Galactic gas-to-dust ratio (e.g., 
iGiiver & Ozelll20^ . and thus the extra dust for the redden¬ 
ing cannot explain the observed X-ray weakness which re¬ 
quires Nu ~ lO^'* cm“^ (Section [4. 111 . This simple testing of 
reddening cannot, of course, exclude the possibility of an in¬ 
trinsically redder spectrum for the X-ray weak WLQs. 

5 . 5 . Spectral Diagnostics of X-ray Weak Quasars 

Studies of the composite SDSS spectra revealed distinct 
characteristics of our X-ray weak PHL 1811 analogs and 
WLQs with respect to the X-ray normal population, including 
redder continua and stronger UV Le II emission. We thus per¬ 
formed a statistical analysis of various emission-line and con¬ 
tinuum properties, searching for possible spectral diagnostics 
of X-ray weak quasars that would help reveal their nature. 

We ran Peto-Prentice tests to assess whether the distribu¬ 
tions of the optical-UV spectral properties of the X-ray weak 
objects differ from those of the X-ray normal objects. We 
combined PHL 1811 analogs and WLQs in these tests con¬ 
sidering that they can be unified (Wll; see also Section 6 
below). We examined the C IV REW, C IV blueshift, C IV 
LWHM, REWs of the A1900, Le II, Le III, and Mg II emis¬ 
sion features, and also the relative SDSS color A{g-i). The 
tests were performed for both the full sample and the C IV 

Recently, IMeusinger & Balafkanl <201411 found that their sample of 
SDSS quasars with weak emission lines have generally bluer continua than 
typical quasars. However, their sample was selected differently from our 
WLQs, and their me thod to generate the composite spectrum also differs from 
ours. Therefore, the lMeusinger & BalafkanI <201411 results cannot be directly 
compared to our findings here. 


The PHL 1811 spectrum is bluer than the SDSS composite spectram 
at short wavelengths (< 1300 A: ILeighlv et alir2007a<l where most of the 
PHL 1811 analogs do not have spectral coverage. 


The SDSS spectrum of the one X-ray normal PHL 1811 analog 
(11537+2716) is more similar to the spectra of X-ray normal WLQs; it is 
bluer and has weaker UV Fe II emission than the composite spectrum of the 


’ http ;//classic.sdss. org/drl / instruments/ imager/index XitaylvIteklWLQs (e.g., see Figure [T^ below). 
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Figure 10. Comparison of the composite SDSS spectra for the (a) 17 X-ray weak PHL 1811 analogs and PHL 1811 and (b) 14 X-ray w eak and 17 X-ray 
normal WLQs, with the y-axis being the flux in arbitrary linear units. In both panels, the SDSS composite spectrum in IVanden Berk et'S]4200lh is also shown for 
compaiison. PHL 1811 and its analogs have similar redder continua than the SDSS composite spectrum, with PHL 1811 having stronger UV Fe II (2250-2650 A) 
emission. The X-ray weak WLQs have generally redder continua than the SDSS composite spectrum, while the X-ray normal WLQ composite spectrum is bluer 
and also exhibits weaker UV Fe II emission. 
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Figure 11. Comparison of the composite SDSS spectra for the (a) 17 X-ray weak PHL 1811 analogs and 8 X-ray weak WLQs and (b) 8 X-ray weak and 10 
X-ray normal WLQs in the C IV subsample. The composite spectra for the X-ray weak PHL 1811 analogs and WLQs are similar to each other. The X-ray weak 
WLQs have redder continua and stronger UV Fe II (2250-2650 A) emission than the X-ray normal WLQs. 


subsample; the results are shown in Table 6. 

We found that statistically, the X-ray weak PHL 1811 
analogs and WLQs have larger Fe II REWs and redder A(g-i) 
colors than the X-ray normal objects, both at the more than 3a 
significance level (3.8 ct and 4.6 ct). The significance levels of 
the differences drop slightly (to 2.9cr and 3.7cr) for the C IV 
subsample, probably due to the smaller sample size. For the 
other spectral properties in either the full sample or the C IV 
subsample, discrepancies are found at the 1.2-2.7cr level: the 
X-ray weak population has in general stronger C IV blueshifts 
(L7cr), larger C IV FWHMs (2.2a), and higher UV emission¬ 
line REWs (e.g., C IV REW at 2.3cr) than the X-ray normal 
population. 

We consider the UV Fe II REW and A(g - i) color as the 
most significant diagnostics of X-ray weak PHL 1811 analogs 
and WLQs. The distributions of the Fe II REWs for the X-ray 
weak and X-ray normal populations are shown in Figure [T3] 
with the results from the Peto-Prentice test listed. The Fe II 
REW measured from the SDSS quasar composite spectrum is 
also plotted for comparison. The X-ray weak objects have in 


fact comparable UV Fe II emission to typical SDSS quasars, 
while the X-ray normal objects have weaker than average Fe II 
emission. This result is co nsist ent with the findings from the 
composite spectra (Section lST^ . We also compared the distri¬ 
butions of the Fe II to C IV REW ratio for the 24 X-ray weak 
and 11 X-ray normal objects with both Fe II REW and C IV 
REW constraints (Figure [T4l i. and no significant difference 
was found (Peto-Prentice test Pnuii = 0.5). 

The distributions of the A(g - i) colors are shown in Fig¬ 
ure with the results from the Peto-Prentice test and 
Kolmogorov-Smirnov (K-S; applicable when there are not 
censored data) test listed. The median A(g - i) color for the 
X-ray weak population is 0.17, and it is 0.01 for the X-ray 
normal population. In the C IV subsample, these values are 
0.23 and 0.03, respectively^ These color offsets are mod¬ 
est, and as pointed out in Section 2.4 of W11, the PHL 1811 

These median color valu es di ffer slightly from those we measured from 
the composite spectra (Section[5i4), as the g-i color of the composite spectra 
(converted to the observed frame using the median redshift) differs from the 
g — i colors of individual objects with different redshifts. 
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Figure 12. Composite SDSS spectrum for the eight X-ray weak WLQs de- 
reddened with Ay = 0.05 mag compared to the composite spectrum for the 
10 X-ray normal WLQs in the C IV subsample. The bottom panel shows the 
ratio of the two spectra. The de-reddening was performed assuming an SMC 
extinction law. The de-reddened X-ray weak composite spectrum matches 
well with the X-ray normal composite spectrum, except for some excess UV 
Fe II and Mg II emission. 

analogs and WLQs are still within th e inclusion area for t he 
SDSS color selection of quasars (e.g.. lRichards et aT]l2002h . 

Figure [16] displays how the X-ray weak and X-ray normal 
sample objects are separated in the relative color vs. Fe II 
REW plot. The normal X-ray emission from J1537-I-2716, 
the only X-ray normal PHL 1811 analog, can be explained by 
its small A(g — i) and Fe II REW (estimated to be 13 ±2 A 
based on the fractional coverage; Table 2) values. There is 
also one obvious outlier in Eigure[T6| J0844H-1245, that has an 
unusually red color, A(g-i) = 0.37, yet is X-ray normal. This 
source has a relatively high radio loudness parameter and it is 
possible that its X-ray emission has some con tribu tion from 
jet-linked emission (see more details in Section l43T l. 

The X-ray weakness of our sample objects is measured by 
the Adox parameter. Due to the significant fraction of sources 
not detected in the X-rays, it is not feasible to perform cor¬ 
relation analysis between Aaox and Ee II REW or A{g-i). 
Instead, we utilized the Kendall’s t and Spearman’s rank- 
order tests in the ASURV package to check whether such a 
correlation exists. These tests were performed on the C IV 
subsample, and the distributions of the quantities are plotted 
in Eigure (TT) The resulting small null-hypothesis probabil¬ 
ities (0.2%-l%) suggest a likely correlation between Aaox 
and Ee II REW or relative color, and the overall trend is such 
that a larger Fe II REW or A(g - i) corresponds to a smaller 
Aaox- We also show the stacked data points in Eigure [17] 
for the 15 Chandra undetected sources (excluding the two 
XMM-Newton undetected sources), which have a mean Aaox 
of -0.99 following the stacking procedure in Section 15.11 
Deeper X-ray observations, converting the Aaox upper limits 
into detections (factors of ~ 10 increase in the exposure times 
are needed given the stacked X-ray flux level), or a larger sam¬ 
ple is required to quantify the possible correlations. 

It is also of interest to probe the correlation between the 
spectral properties and aox, which is an indicator of the 
X-ray flux level. In the C IV subsample, the objects have 
a narrow distribution in L2.S00 A (within about an order of 
magnitude) and thus a na rrow range of ex pected aox values 
(within ^0.14) from the lJust et all (l2007h aox-L^^np ^ rela- 
tion. Therefore, replacing Aaox with aox in Eigure llTl does 


not change the distributions significantly, and similar trends 
still exist in the sense that a larger Ee II REW or A{g - i) cor¬ 
responds to a smaller aox- 

5.6. Selection Bias in the X-ray Weak Quasar Diagnostics? 

We found that X-ray weak PHL 1811 analogs and WLQs 
generally have larger Fe II REWs, as well as larger C IV 
blueshifts, C IV EWHMs, and UV line (C IV, the A1900 com¬ 
plex, Ee III, and Mg ll) REWs at a less significant level, than 
the X-ray normal population. However, PHL 1811 analogs 
and WLQs were mixed to form a large sample for statistical 
analysis, which could in fact introduce a selection effect due 
to the different selection criteria for these two types of objects. 
As detailed in Section 2, the PHL 1811 analogs were selected 
with a less stringent requirement on the C IV REW, and no re¬ 
quirements on the REWs of the A1900 complex and Mg II, but 
with extra requirements on the C IV blueshift and UV Ee II and 
Ee III strength. The fraction of X-ray weak quasars is signif¬ 
icantly higher among the PHL 1811 analogs than the WLQs, 
and X-ray weak PHL 1811 analogs also constitute a substan¬ 
tial fraction (> 50%) of the X-ray weak population analyzed 
here. A consequence is that all the above criteria would stand 
out as spectral diagnostics of X-ray weak quasars in our anal¬ 
ysis while only one or some of these might really be respon¬ 
sible for the weak X-ray emission. 

To avoi d this possible selection bias, we performed the tests 
in Section [ 53 ] using only the WLQs in the C IV subsample, 
which has eight X-ray weak and 10 X-ray normal objects. Not 
all of these objects have measurements of all the emission¬ 
line properties. Only three properties in Table 6 differ at the 
more than 1.5cr significance level between the X-ray weak and 
X-ray normal populations: Ee II REW (2.0cr), Mg II REW 
(L6 ct), and A(g-i) {2.1a). The less-significant test results 
could either be due to the selection bias being removed or 
simply the smaller sample size. The latter effect is evident 
in the A(g-i) results {3.1a to 2.1a), which are not affected 
by the selection bias. Nevertheless, Ee II REW and A{g-i) 
remain the most robust diagnostics of X-ray weak quasars. 

It is also possible that the relatively larger REWs of the UV 
lines (e.g., C iv) of the X-ray weak population are real instead 
of being a selection effect, if they are intrinsically connected 
with the relatively larger Ee II REW, e.g., in the shielding-gas 
scenario discussed in Section 1631 below. With our currently 
limited sample of X-ray observed WLQs, it is difficult to as¬ 
sess whether there is a selection bias in the diagnostic analy¬ 
sis. 

6. DISCUSSION 

6.1. Unification of PHL 1811 Analogs and WLQs with the 
Shielding-Gas Scenario 

The overall similarities between the X-ray weak PHL 1811 
analogs and WLQs, including their IR-UV SEDs (Figures [7] 
and[3i, composite SDSS spectra (Figure [TT]i, and degrees of 
X-ray weakness (Figure [5]l, suggest that PHL 1811 analogs 
and WLQs are physically similar types of object. PHL 1811 
analogs are empirically a subset of WLQs (with small dif¬ 
ferences causedby the differing selection approaches; Sec¬ 
tions O and O, and with the additional requirements of 
PHL 1811-like emission-line properties, we preferentially se¬ 
lected X-ray weak WLQs as the PHL 1811 analogs. Give n our 
analysis of the X-ray weak quasar diagnostics in Section 1531 
the criterion of strong UV Fe II emission most likely influ¬ 
enced this outcome, while the criteria of large C IV blueshift 
and strong UV Fe III emission might also play a role. 
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Fe II REW (A) 



Fe II REW (A) 


Figure 13. Distributions of the Fe II REWs for the (a) full sample and (b) C IV subsample. The hatched red and solid blue shaded histograms and arrows 
(representing limits) are for the X-ray weak and X-ray normal objects, respectively. We list the probability of the two distributions being drawn from the same 
p arent population based on the Peto-Prentice test. We also list, and plot as the vertical dashed line, the Fe II REW measured from the SDSS composite spectrum 
in IVanden Berk et all <200111 . The X-ray weak objects have on average larger Fe II REWs than the X-ray normal objects. 
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Fe II REW/C IV REW 

Figure 14. Distribution of the Fe II REW to C IV REW ratios for the X-ray 
weak (red) and X-ray normal (blue) objects. The X-ray weak and normal 
objects have similar distributions of these line ratios, and the Peto-Prentice 
test also suggests no significant difference (Pnull = 0.5). 

Based on the properties of the small pilot sample of 
PHL 1811 analogs, W11 proposed a simple shielding-gas sce¬ 
nario to unify PHL 1811 analogs (X-ray weak WLQs) and 
X-ray normal WLQs. In this model, the shielding gas has a 
sufficiently high covering factor to shield all or most of the 
BELR from the ionizing continuum, resulting in the observed 
weak UV emission lines. Whether we observe an X-ray weak 
or X-ray normal quasar is then an orientation effect, depend¬ 
ing upon whether our line of sight intersects the X-ray absorb¬ 
ing shielding gas. X-ray normal WLQs are observed at small 
inclination angles, while X-ray weak WLQs and PHL 1811 
analogs are observed at larger inclination angles. 

With our extended sample of PHL 1811 analogs and WLQs, 
the spectral stacking results suggest that the X-ray weak ob¬ 
jects are generally ab sorb ed instead of being intrinsically 
X-ray weak (Section EB. The presence of the signifi¬ 
cant fraction of X-ray normal WLQs (Ri 50%) also indi¬ 
cates that the weak UV line emission in PHL 1811 analogs 
and WLQs cannot universally be attributed to intrinsic X-ray 
weakness. The average X-ray weakness factor for our X-ray 
weak PHL 1811 analogs or WLQs in the C IV subsample is 


r; 40, corresponding to Vh « 9 x 10^^ cm ^ (Section 14.1b . 
The spectral analysis of J1521-F5202 also revealed strong 
X-r ay ab sorption with Nn = (1.26;|;o 4 g) x 10^^ cm“^ (Sec¬ 
tion [Q l. We note that the actual Nh values could be much 
larger if the observed X-ray emission is dominated by re- 
flected/scattered X-rays r ather than direct transmission (e.g., 
iMurphv & Yaaoobll2009l) : a substantial contribution from re¬ 
flected/scattered X-rays is expected for such levels of X-ray 
weakness based on observations of, e.g., local Seyfert 2 galax¬ 
ies. These findings provide further support for the Wll 
shielding-gas scenario, where a soft ionizing continuum due 
to small-scale absorption is the key for creati^ the BELR line 
properties of PHL 1811 analogs and WLQsEj 

It has been suggested that the EWHMs of low-ionization 
lines (e.g., H/3, Mg ll) have an orie ntation dependence due 
to a flattened BELR geometry (e.g., IWills & Browne! 119861 : 
[Runnoe et alJl201^IShen & Hdl201^ , with the EWHMs gen¬ 
erally being larger at larger inclination angles (although there 
is considerable object-to-object scatter). There are 23 X-ray 
weak objects and 17 X-ray normal obj ects in our sample hav¬ 
ing Mg II EWHM measurements in the lShen et alJ (1201 Ih cat- 
alogl3 The X-ray weak objects have relatively larger Mg II 
EWHMs than the X-ray normal objects (mean EWHM values 
5400 ± 600 km s“* vs. 4100 ± 600 km s“' and median values 
5300 km s“^ vs. 2900 km s“Q, consistent with the shielding- 
gas scenario where the X-ray weak objects are viewed at 
larger inclination angles. However, a K-S test between the 
two sets of Mg II EWHMs does not indicate a highly signif¬ 
icant difference, so we consider this result only suggestive at 
present. 

One possible consequence of the shielding-gas scenario is 
that the Aaox distribution of PHL 1811 analogs and WLQs 
would appear bimodal, at least to hrst approximation, as the 

The He II AI640 line has been suggested to be a good indicator of 
the strength of the ionizing c ontinuum, with a smaller REW implying a 
weaker/softer continuum (e.g., ILeighlvlf2Q04l: IBaskin et alJI2Ql^. We mea- 
sured the He II REWs for our sample objects the same way as in iBaskin et alj 
(Ml), and the results do show weaker than average He II REWs for our 
PHL 1811 analogs and WLQs. 

Upon visual inspection, we consider that it is necessary to subtract a 
narrow Mg II line component for J1629+2532, and the EWHM value of the 
broad component was adopted. For the other objects, the EWHM values of 
the whole Mg II profile are adopted. 
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Figure 15. Distributions of the relative color, A(g —i), for the (a) full sample and (b) C IV subsample. The hatched red and solid blue shaded histograms 
represent X-ray weak and X-ray normal objects, respectively. The vertical dashed line indicates a relative color of zero. The probabilities of the two distributions 
being drawn from the same parent population based on the Peto-Prentice test and K-S test are listed. The X-ray weak objects are in general redder than the X-ray 
normal objects and typical SDSS quasars. 
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Figure 16. Relative color vs. Fe II REW for the (a) full sample and (b) C IV subsample, including the X-ray weak PHL 1811 analogs (magenta squares). X-ray 
normal PHL 1811 analog (green diamond). X-ray weak WLQs (red circles), and X-ray normal WLQs (blue triangles). The X-ray normal PHL 1811 analog has 
an estimated Fe II REW of 13 ± 2 A based on the fractional coverage of Fe II (Table 2). The dashed lines represent the median SDSS color an d the Fe II REW in 
the SDSS composite spectrum. We also mark two objects, J0844-I-1245 and J1156+1848, that ar e con sidered radio intermediate (R > 10) in the lShen et al.l 1201 II) 
SDSS quasar catalog using a slightly different definition of R from the one used here (Section [4.3i . The WLQs and PHL 1811 analogs that have redder colors 
and/or larger Fe II REWs than the average SDSS values are more likely to be X-ray weak quasars. 


sources are either heavily X-ray absorbed or X-ray normal. 
Due to the large fraction of undetected sources in the X-ray 
weak population, the Aaox distributions for our samples 
(Figure |6]l cannot reveal clear bimodality. However, the very 
large average X-ray weakness factors (> 100; Section 14. Il l 
for these undetected sources do suggest that the Aaox values 
cannot have a continuous distribution and bimodality is plau¬ 
sible. Deeper X-ray observations are required to constrain 
better the Aaox distribution. 

Based on the X-ray and multiwavelength properties of 
PHL 1811 analogs and WLQs, some basic physical require¬ 
ments on the shielding gas can be obtained; 

1. The shielding gas has a large column density of X-ray 
absorption (at least 10^^ cm“^ and likely much larger). 

2. It lacks accompanying C TV BALs or mini-BALs at 
least along our line of sight. 

3. It lies closer to the SMBH than the BELR, so that it 


can screen the nuclear EUV and X-ray emission from 
reaching the BELR. Eurthermore, it should have a large 
covering factor to the BELR. 

4. The “waste heat” from the absorbed high-energy emis¬ 
sion (that otherwise would have reached the BELR) is 
likely re-emitted in the unobserved EUV as the contin¬ 
uum IR-UV SED appears normal. Such re-emission 
in the EUV is indicative of a small-scale absorber on a 
scale of m lORs (Rs = is the Schwarzschild 

radius). 

Given these properties, the shielding gas might naturally be 
understood as a geomet rical ly thick inner accretion disk (e.g., 
a slim disk; see Section lh^ belowl. This would require much 
of the BELR gas to be in an equatorial configu ration, as sup¬ 
ported by observations (e.g., IShen & Holl2014l) . 

6 .2. A Geometrically Thick Disk Scenario for the X-ray 
Absorbing Shielding Gas 














































16 


0.4 

0.2 

0.0 

-0.2 


o 

< -0.4 


- 0.6 

- 0.8 

- 1.0 


(o) 


o 

o o 


C IV subsomple 

Kendoll's T test: 

Pn„„ = 0.3% 

Speormon's ronk-order test: ’ 
Pn„„ = 0.4% 




b 

1 iU° 1 i 


stocked mean of 15 upper limits 


20 40 

Fe II REW (A) 


60 



0.4 

(b) 

-^^-(- 

1 

-n -^^ ^^ - 1-,-^-,- 




1 

Q C IV subsomple 


0.2 

- 

O 1 

° Kendoll's t test: 




01 

P„„„=0-2% 


0.0 

- 

'o 

° Speormon's ronk-order test 




I (55 

o ' 



-0.2 


”1 “ 

0 

i 

o 

0 

< 

-0.4 

- 

1 

1 

1 

•1 ^ i i O o °| 1 

4 i 




f i 

i 0 0 i 


-0.6 

- 

1 

o 


-0.8 

- 







Stocked meon of 15 upper limits 


-1.0 

- 


★ 


-0.2 0.0 0.2 0.4 

Relative Color A(g-i) 


0.6 


Figure 17. Aaox vs. (a) Fe II REW and (b) relative color for the C IV subsample. The horizontal gray dashed lines represent the division between X-ray weak 
and X-ray normal quasars in this study. The vertical gray dashed lines display the Fe II REW in the SDSS composite spectrum and the median SDSS color, 
respectively. Eor the undetected 15 objects with Chandra data (the other two undetected sources have XMM-Newton observations) in this figure, we show their 
stacked Aoqx value and the mean Ee II REW (or relative color) as the magenta star. We performed both the Kendall’s r and Spearman’s rank-order tests in 
the ASURV package for correlation analyses. These tests work with data having upper and/or lower limits, and the resulting small null-hypothesis probabilities 
(0.2%-!%) suggest significant coiTelations between Aaox and Fe II REW or relative color. 


Although the X-ray weak and X-ray normal PHL 1811 
analogs and WLQs can be unified under the W11 shielding- 
gas scenario, the physical nature of this shielding gas, that can 
produce such quasars with extreme emission-line and X-ray 
properties, remains mysterious. Base d on its physical re¬ 
quirements enumerated above (Section fe.lb . we propose that 
a geometrically thick inner accretion disk may naturally serve 
as the shielding gas. We discuss below such a scenario in 
the context of a puffed-up disk due to rapid or even super- 
Eddington accretion. However, we first note that the basic 
idea of our model stands irrespective of the exact physical 
processes leading to the geometrical thickness of the disk. 

Given the small fraction of PHL 1811 analogs and WLQs 
among SDSS quasars, a geometrically thick accretion disk 
with a sufficiently large scale height to shield the BELR al¬ 
most fully should be a rare phenomenon. Under our selec¬ 
tion criteria for the PHL 1811 analogs, 66 candidates were 
identified from Ri 7200 SDSS quasars satisfying the redshift, 
mag nitu de, radio-loudness, and non-BAL requirements (Sec¬ 
tion EB- The corresponding fraction is r; 0.9%, consistent 
with the estimation of ;< 1.2% in W11. The fraction of WLQs 
among SDSS quasars is likely larger by a factor of R 2 due 
to the addition of the X-ray normal population. The rarity of 
these quasars suggests a link to some extreme physical prop¬ 
erty. One relevant physical quantity is the Eddington rati o, 
as suggested by earl ier studies (e.g.. iLeighlv et aHl2007allbt 
iShernmer et al.i2010l) . 

PHL 1811 itself has an estimated Eddington ratio of 
R 1.6 (ILeighlv et al.l l2007bl) . PHL 1811, along with sev¬ 
eral PHL 1811 analogs and WLQs with rest-frame opti¬ 
cal spectra (J15214-52 02, 2QZ J2154-3056 , and the two 
high-redshift WLQs in IShernmer et al.ll2010l) . have weak or 
undetected [O III] A5007 narrow emission lines, sugges¬ 
tive of high Eddi ngton ratios (e.g.. iBoroson & Greenlll9^ 
IShen & Hal2014l) . Moreover, several studies have found that 
as LBoi/TEdd increases, the C IV REW genera lly decreases 


and 

the C IV blueshift also increases (e.g.. iBachev et al. 

2004 

; [Baskin & Laoi1l2004l [Richards et al.ll201 iT Shen & Ho 

2014 

ilSulentic et al.l2014HShemmer & Liebeit20L5|). In these 


correlations, there is only limited sampling in the super- 
Eddington or low (< 10 A) C IV REW regime, but the overall 


trends suggest that the Eddington ratio grows as one moves 
from typical quasars toward WLQs in the C IV REW vs. 
blueshift space (Eigure[^). By the nature of our selection of 
the PHL 1811 analogs and WLQs, demanding that the BELR 
does not produce normal emission lines, we may have recov¬ 
ered effectively a population of quasars with (extremely) high 
Eddington ratios. It is difficult to measure LBoi/TEdd directly 
for our exceptional objects, as the SMBH masses estimated 
from the line-based virial method are likely high ly uncertain 
and perhaps systematically in error (Section EB- However, 
based on the empirical T-Laoi/TEdd relations, our joint spec¬ 
tral anal ysis of the X-ray normal subsample in Sections 15.21 
and 15.31 does indicate a high Eddington ratio (Luoi/TEdd ~ 1) 
in general for our quasars 

A high Eddington ratio is naturally connected to our 
proposed geometrically thick accretion disk. When 
the Eddington ratio is high (LBoi/Tedd ^ 0.3), optically 
thick advection becomes important and a slim accre¬ 
tion d isk (e.g.. I A bramo wicz et al.l 119881: [Mineshig e et alJ 


2000t IWang & Netzed 120031: lOhsuga & Mineshigel 1201 1 : 
Straub et al.ll201 U iNetzer & Trakhtenbroj 120141 : IWang et alJ 


2014ai and reference s therein) is a more appro priate solu 


tion than the standard fShakura & SunvaeB (119731) thin accre¬ 
tion disk. The slim disk has a geometrically thick inner re¬ 
gion which, in the case of PHL 1811 analogs and WLQs, 
might act as the shielding gas that blocks the nuclear ioniz¬ 
ing continuum from reaching the BELR (e.g., lLeighlvl[2004l : 
IWang et ani2014bh . Recent three-dimensional global MHD 
simulations of super-Eddington accre tion disks with ac curate 
and self-consistent radiative transfer (l.liang et alJl2014l : Y.-E. 
Jiang 2015, private communication) or MHD simulations of 
super-Eddington accretion d isks accounting for ge neral rela¬ 
tivity and magnetic pressure (ISadowski et al.ll^l4^ also show 


Narrow-line Seyfert I galaxies (NLSls) also generally have steep X-ray 
spectra and high Ed dington ratios. PHL 1811 itself is considered a NLSl 
jLeighlv et al.ir2QQl1) . and NIR spectroscopy of a limited sample of WLQs 
shows that they have in general narrow and strong opticd Fe II emis¬ 
sion (see Section 5.1 of IPlotkin et alj[2Q15D . similar to NLSls. However, 
as the line widths and REWs have likely dependences on luminosity, our 
PHL 1811 analogs and WLQs are not directly comparable to local, less- 
luminous NLSls. 
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this geometrically and optically thick inner region. An ex¬ 
tremely high Eddington ratio is probably required for the in¬ 
ner disk to be puffed up sufficiently to reach the large neces¬ 
sary covering factor to the BEER. This thick inner disk might 
also block a significant portion of the ionizing radiation from 
reaching the narrow li ne region, leading to the observed weak 
[O III] emission (e.g.. [Boroson & Greenlll9^ . A schematic 
illustration of this scenario is shown in Eigure (TS] In or¬ 
der for the puffed-up disk to block the nuclear high-energy 
emission from reaching the BEER, the disk corona must be 
fairly compact. Studies of the rapid X -ray yariability (e.g., 
IShemmer et al.l20' ^lUttley et al.l201 ^ and X -ray microlens- 
ing (e.g. jDai et al~ 120 lOt [Morgan et al.ll2012h of AGNs haye 
constrained the sizes of the X-ray emitting regions to Ri SR^ in 
these systems, consistent with the requirements of our model. 

The expected SED from a super-Eddington accretion disk 
is uncertain. In the optical-UV, it likely has a power-law 
shape similar to a standard thin accretion disk, but in the 
EUV to s oft X-ray range (r! 10 -100 eV), it is probably flat¬ 
ter (e.g.. IWang & Netzerl 120031) . The nuclear high-energy 
emission absorbed by the geometrically thick inner disk (the 
shielding gas) will likely enhance its EUV (e.g., r; 50 eV) 
emission, and thus the source will appear to haye a typi¬ 
cal quasar continuum SED from the IR to UV (Section [4.21 l. 
The column density through the puffed-up inner accretion 
disk is high (Nh ^ 10^"^ cm“^), much larger than the ayer- 
age Nh of Ri 9 x 10^^ cm~^ roughly estimated for our X-ray 
weak objects (Section 16. 11 1. Therefore, the obseryed X-ray 
emission from these objects is likely dominated by indi¬ 
rect, reflected/scattered X-rays instead of transmitted emis¬ 
sion through the disk. The nuclear X-ray emission could be 
Compton reflected by the accretion disk and/or scattered by 
an ionized scattering medium on a larger scale than the in¬ 
ner disk. A broad range of X-ray weakness factors (e.g., 
Eigure |6ll might be expected as the reflection/scattering effi¬ 
ciency yaries for indiyidual objects, but the oyerall leyels of 
X-ray weakness should be high for reflection/scattering dom¬ 
inated spectra so that bimodality in t he Aaox distribution is 
still likely expected (see Section lOT i. If the obseryed X-ray 
emission is dominated by reflection, a strong Ee Ka emission 
line with a REW of order 1-2 keV is usually expected (e.g., 
iGhisellini et aklflO^lMatt et al.ll 19961) . which cannot be con¬ 
strained with our limited data. 

If our line of sight is close to the edge of the inner bulge 
of the puffed-up disk in some objects, a small change in 
the coyering factor of the inner disk due to, e.g., accretion- 
disk instability could result in transitions between X-ray 
weak and X-ray normal states on timescales of years. Such 
a phenomenon may haye been obseryed in the z = 0.396 
quasar PH L 1092. It has simila r emission-line properties to 
PHL 1811 (iLeighly et al.ll200^ W11), and it yaried between 
an X-ray normal quasar and an X-ray weak quasar oyer a 
timescale of years with a maximum yariability factor of r; 260 
in its 2 keV flux; meanwhile, its UV emission lines did not 
show such a drastic change (iMiniutti et al.l l2009l 1201 2^ as 
would be expected if the coyering-factor change is small. Eur- 
ther X-ray monitoring obseryations of PHL 1811 analogs and 
WLQs are required to constrain the frequency and duration of 
such transitions, which may be a useful probe of the geomet¬ 
rically thick disk scenario. 



Figure 18. A schematic diagram of the geometrically thick disk scenario for 
PHL 1811 analogs and WLQs. When the accretion rate is very high, the in¬ 
ner region of the accretion disk is puffed up significantly, which serves as the 
shielding gas in the W11 scenario. It blocks the nuclear ionizing continuum 
from reaching the BELR, resulting in the observed weak line emission. The 
dashed lines illustrate a standard thin accretion disk, where the BELR is illu¬ 
minated normally. The inner region of the disk can also absorb the coronal 
X-ray emission, resulting in an X-ray weak PHL 1811 analogAVLQ when 
the disk inclination angle is moderate or an X-ray normal object when the 
inclination angle is small. A standard accretion disk wind is launched from 
the outer region of the disk. The BELR contains b oth a static componen t and 
a wind component, at least for the C IV line (e.g.. lRichards et ^1201 Ih . We 
annotate the approximate scale of the disk in several places, in units of the 
Schwarzschild radius (Rs). 


Based on our geometrically thick disk scenario for 
PHL 1811 analogs and WLQs, these extreme and rare quasars 
are likely not a distinct population, but are instead extreme 
members of the continuous population of quasars. The shield¬ 
ing effect from a puffed-up inner disk likely exists beyond 
these extreme objects, and it might be applicable at a milder 
level to quasars with lo wer Edding t on rat ios; in either the 


_ agt on rat ios 

slim-disk model or the iJiang et ^ (l2014h simulations, the 
disk is unlikely as thin and flat as what a standard disk model 
describes as long as LBoi/TEdd ^ 0.3. When the Eddington ra¬ 
tio is smaller than those of PHL 1811 analogs and WLQs, the 
radius and scale height of the puffed-up disk decreases, and 
its covering factor to the BELR also decreases, leading to a 
larger C IV REW than those of PHL 1811 analogs and WLQs. 

Given a continuously varying covering factor of the in¬ 
ner disk that is LBoi/TEdd dependent, a substantial fraction 
of the quasar C IV REW distribution (spanning about two 
orders of magnitude; Eigure [^) could be governed by the 
Eddington ratio. Additional factors shaping this distribu¬ 
tion include anisotropic continuum emission, anisotropic line 
emission, gas metallicity, C IV BELR geometry (setting the 
fraction of ionizing radiation that is ca ptured), and poten 


tial self-shadowing of the BEER (e.g., [Koristae^^ 


Korista & (joadll2004t iGoad et akl 120 121 [Baskin et ~ 


1998; 


2014; 


IWang et alJl2014bh . Eurthermore. shielding of the nuclear X 
ray emission can prevent the disk wind from b eing overion¬ 
ized and help the launc hing of the wind (e.g., iMurrav et alJ 
119951 IWoga et al.ll2000li . A larger covering factor of the geo¬ 
metrically thick disk (larger TBoi/TEdd) could thus produce a 
stronger outflowing wind and result in a larger C IV blueshift. 
The distribution of quasar C IV blueshifts (Eigure [^) could 
thus also be partly explained by the disk shielding scenario. 


6.3. Broader Implications of Geometrically Thick Disks for 
Quasar Emission Lines 


6.4. X-ray Weak Quasar Diagnostics and the Shielding-Gas 

Scenario 
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We investiga ted t he X-ray weak quasar diagnostics discov¬ 
ered in Section l531 and found that they can be plausibly fitted 
into the W11 shielding-gas scenario, where the X-ray weak 
objects are viewed at larger inclination angles than the X-ray 
normal objects. 

The redder A(g - /) color of the X-ray weak population 
can be explained by m ild excess intrinsic reddening {Ay ~ 
0.05 mag; Sectionif the dust tends to reside in the equa¬ 
torial plane of the quas ar (e.g., in a dust-driven outer wind 
scenario as discussed bv lElvisll2012h which leads to more ex¬ 
tinction at larger inclination angles. Another possible scenario 
is that the accretion-disk emissio n is n ot isotropic (e.g., from 
a puffed-up slim disk; see Section l672l) . and less UV emission 
is received at larger inclination angles, producing an intrinsi¬ 
cally redder continuum. 

The larger REWs of the UV lines (C IV, the A1900 com¬ 
plex, Ee III, Ee II, and Mg ll) in the X-ray weak population 
compared to those in the X-ray normal population may be a 
consequence of the aspect dep e ndent effects of accretion-disk 
emission (e-g-, iNetzerl 119901: iNetzer & Trakhtenbrolj 12014t 
IWang et al.l|2014bl) . The continuum flux from a standard ac¬ 
cretion disk scales approximately as cos i(l+2 cos /), where 
i is th e inclination angle of the disk (iNetzer & Trakhtenbrod 
120141) . In the case of a slim disk, the disk continu um emission 
is sign ificantly anisotropic (e.g., see Eigure 4 of IWang et alJ 
12014bl) . The emission-line flux, on the other hand, is presum¬ 
ably from a BEER structure which overall emits more isotrop¬ 
ically. The resulting line REW thus depends on the inclination 
angle and may be increased by a factor of a few moving from 
a small angle for an X-ray normal WLQ to a moderate angle 
for an X-ray weak WLQ. In this scenario, the ratios of the 
emission-line REWs for the X-ray weak and X-ray normal 
populations should not differ significantly. This interpretation 
is supported by the distributions of the Ee II to C IV REW ra¬ 
tios for the two populations in our sample (Eigure O, which 
show no significant difference (Peto-Prentice test Pnuii = 0.5). 

The theoretical angular dependence of disk continuum 
emission is slightly stro nger toward l onger wavelength in the 
UV (e.g., Eigure 32 of lNetzeJll990l) . The Ee II emission is 
at a longer wavelength than the other studied UV lines except 
Mg II, and it is also relatively stronger than the other weak 
UV lines in WLQs (leading to more reliable measurements). 
These two factors combined probably cause the Ee II REW 
to be the best X-ray weak quasar diagnostic among all the 
REWs in our analysis (C IV, the A1900 complex, Ee III, Ee II, 
and Mg ll). 

The other two less-significant diagnostics of X-ray weak 
quasars, the C IV blueshift (1.7 ct) and EWHM (2.2cr), can be 
generally incorporated in the W11 scenario considering that 
the C IV line has an accelerating disk-wind origin so that the 
line may have both a larger blueshift and a larger velocity 
broadening if our line of sight is more aligned with directions 
along which the wind is strongly accelerated. A large C IV 
blueshift was one of the selection criteria for our PHL 1811 
analogs, but it appears less efficient than the requirement of 
large UV Ee II REW for selecting X-ray weak quasars; the 
orientation-angle dependence of the line-of-sight wind veloc¬ 
ity is probably less significant than the orientation-angle de¬ 
pendence of the line REWs. 


7. SUMMARY AND FUTURE WORK 
7.1. Summary of Main Results 


In this paper, we have presented X-ray and multiwavelength 
investigations of the nature of a large sample of PHL 1811 
analogs and WLQs. The key points from this work are the 
following: 

1. We obtained Chandra exploratory (1.5-9.5 ks) observa¬ 
tions of 10 PHL 1811 analogs and 22 WLQs. We also 
acquired a 40 ks Chandra observation of J1521H-5202, 
a PHL 1811 analog that is one of the most luminous 
SDSS quasars. We measured their X-ray photomet¬ 
ric properties, and performed basic spectral analysis for 
J1521H-5202 which suggests strong intrinsic X-ray ab¬ 
sorption. Including the previous samples in W11 and 
W12, we constructed a large sample of 18 PHL 1811 
analogs and 33 WLQs at z = 0.5-2.9 that have X-ray 
observations. See Sections 2 and 3. 

2. Out of the 18 PHL 1811 analogs, 17 (94%) are X-ray 
weak; out of the 33 WLQs, 16 (48%) are X-ray weak. 
The selection criteria of PHL 1811-like emission-line 
properties worked effectively for finding X-ray weak 
quasars, and we can accommodate the one X-ray nor¬ 
mal PHL 1811 analog based on its SDSS spectral prop¬ 
erties. The average X-ray weakness factor for the X-ray 
weak PHL 1811 analogs or WLQs in the C IV subsam¬ 
ple is Ri 40. The Aaox distributions of the PHL 1811 
analogs and WLQs are significantly different from that 
of typical quasars. See Section 4.1. 

3. We constructed IR-X-ray continuum SEDs for the 
PHL 1811 analogs and WLQs. Both the X-ray weak 
and X-ray normal groups have IR-UV continuum SEDs 
similar to those of typical quasars. See Section 4.2. 

4. The stacked effective power-law photon indices for the 
X-ray weak subsamples are relatively hard (Eeff ~ 0.8- 
1.5). These, together with the strong X-ray absorp¬ 
tion found in J1521H-5202, suggest that the X-ray weak 
PHL 1811 analogs and WLQs on average are X-ray 
absorbed, although the possibility of intrinsic X-ray 
weakness cannot be completely excluded for some ob¬ 
jects. See Section 5.1. 

5. We performed joint spectral fitting for the 18 X-ray nor¬ 
mal PHL 1811 analogs and WLQs. The best-fit hard 
X-ray photon indices are r = 2.18 ±0.09 for the 18 ob¬ 
jects and r = 2.26 ±0.11 for the 10 objects in the C IV 
subsample, suggesting a high Eddington ratio in gen¬ 
eral for these X-ray normal objects (and also the X-ray 
weak objects as they can be unified). See Sections 5.2 
and 5.3. 

6 . We compared composite SDSS spectra for the X-ray 
weak and X-ray normal PHL 1811 analogs and WLQs, 
and investigated their optical-UV spectral properties as 
diagnostics for identifying X-ray weak quasars. Statis¬ 
tically, the X-ray weak PHL 1811 analogs and WLQs 
have significantly (> 3a) larger UV Ee II REWs and 
redder A(g - i) colors than the X-ray normal popula¬ 
tion. They also have in general larger C IV blueshifts 
and EWHMs, and larger UV line (C IV, the A1900 
complex, Ee III, and Mg ll) REWs at a less signifi¬ 
cant level, than the X-ray normal population. The nor¬ 
mal X-ray emission from J1537H-2716, the only X-ray 
normal PHL 1811 analog, can be understood given its 
small A(g — i) and Ee II REW. See Sections 5.4-5.6. 
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7. The PHL 1811 analogs empirically appear to be a 
subset of WLQs in general, and these two groups 
of objects can be unified under the Wll shielding- 
gas scenario. Due to the additional requirements of 
PHL 1811-like emission-line properties, we preferen¬ 
tially selected X-ray weak WLQs as the PHL 1811 
analogs. The X-ray absorption found in J1521H-5202 
and our stacking results provide further support for the 
Wll shielding-gas scenario. See Section 6.1. 

8 . With the requirement of small C IV REWs for 
PHL 1811 analogs and WLQs, we may have selected 
effectively a population of quasars with geometrically 
thick inner accretion disks that can block the ionizing 
continuum from reaching the BELR and act naturally 
as the shielding gas of the Wll scenario. These quasars 
probably have unusually high Eddington ratios so that 
the inner disk is significantly puffed up (e.g., a slim 
disk). See Section 6.2. 

9. PHL 1811 analogs and WLQs are likely not a distinct 
population, but are extreme members of the continu¬ 
ous population of quasars. Shielding of the BELR by 
a geometrically thick disk is thus perhaps generally ap¬ 
plicable to quasars with lower Eddington ratios, gov¬ 
erning substantial fractions of the observed broad dis¬ 
tributions of quasar C IV REWs and blueshifts shown 
in Figure!^. See Section 6.3. 

7.2. Future Work 

Further investigations are needed to test the ideas in this 
study. We suggested a geometrically thick disk scenario 
to explain the nature of PHL 1811 analogs and WLQs and 
to unify these objects. If these quasars are indeed systems 
with unusually high Eddington ratios, we would expect a 
larger fraction of them at higher redshifts, as the quasar 
Eddington ratio generally grows as redshift incr eases (e.g., 
iNetzer & Trakhtenbrotil200^lShen & K^IMTl) . WLQs do 
appear to be more common at higher redshifts from studies 
of limited WLQ samples at z~ 3-6 (iDiamond-Stanic et alJ 
l2009tlBanados et al.ll2014l) . To test further this notion, a large 
sample of WLQs selected systematically across a broad range 
of redshift is required. 

However, there is no simple way to define universal selec¬ 
tion criteria for WLQs at all redshifts based solely on SDSS 
spectroscopy due to the limited spectral coverage (e.g.. Sec¬ 
tion 3.2 of W12). Broader spectral coverage, such as NIR 
spectra for high-redshift WLQs and UV spectra for low- 
redshift WLQs, are needed to study the correlations of weak 
emission lines and thus define consistent selection criteria 
for WLQs at different redshifts. Currently, only a small 
sample of WLQs have such broad spe ctral coverage (e.g., 
IShemmer et al.ll20T^ iPlotkin et al.ll2oi5T) , and NIR/UV spec¬ 
troscopic observations of more WLQs are needed. Moreover, 
given the unified nature of PHL 1811 analogs and WLQs, the 
WLQ selection criteria (REW < 5 A for all emission features) 
should be relaxed to C IV REW < 10 A, which will allow the 
selection of a WLQ sample that is more consistent with the 
selection of PHL 1811 analogs. 

If the distributions of quasar C IV REWs and blueshifts 
(Figure [^) are indeed partly the result of the varying cov¬ 
ering factor of a geometrically thick disk (which depends on 
the Eddington ratio) as we suggest in Section 6.3, we would 
also expect to see more quasars in general having smaller C IV 


REWs and larger C IV blueshifts at higher redshifts due to 
the redshift evolution of the quasar Eddington ratio. Caution 
should be applied when performing such analyses as the C IV 
REWs and blueshifts might have other luminosity and/or red¬ 
shift dependences (e.g., the Baldwin effect). 

The most- luminous quasy at z > 6 discovered recently is 
also a WLQ (IWu et al.ll2015h . consistent with our expectation 
that the WLQ fraction rises with redshift. It probably has a 
high or even super-Eddington accretion rate, as we propose 
for WLQs (Section 6); a luminosity exceeding the Edding¬ 
ton limit may help explain its very high luminosity. This 
quasar has a M g ll-based virial mass of Ri 1.2 x 10'° Mq 
(IWu et al.ll2015h . However, as we discussed in Section 5.3, 
virial-mass estimates for these extreme quasars are highly un¬ 
certain and perhaps systematically in error (e.g.. Mg ll-based 
virial masses are on average « 3 times larger than H/3-based 
virial masses for five of our objects); this factor could allevi¬ 
ate the challenge of growing such a massive black hole in the 
early universe. A Chandra observation has been scheduled 
for this WLQ (PI: X. Fan). Our X-ray results on WLQs sug¬ 
gest there is a r: 50% chance that this quasar is X-ray weak. 
This chance is likely higher based on this object’s relatively 
red co ntinuum, with a spectral index of ax = -1.43 (IWu et al.l 
1201 5h as compared to the average spec tral index of ax = -1.7 2 
for SDSS quasars (e.g.. Section 3.2 of iKrawczvk et al.ll20L5h . 
Quantification that its UV Fe II emission is strong would in¬ 
crease this chance further (see Section 5.5). 

Due to the substantial uncertainties in the SMBH mass esti¬ 
mates, it is challenging to measure quantitatively the Edding¬ 
ton ratios of PHL 1811 analogs and WLQs. However, the Ed¬ 
dington ratio appears to be the primary driver of quasar Eigen¬ 
vector 1 (e.g..lBc)roson & Greenll99^lNetzer & Trakhtenbrotl 
120071 IShen & Hoi 120141) . which is dominated by the opti¬ 
cal Fe ll/H/3 ratio, [O III] REW, and H/3 FWHM. NIR spec¬ 
troscopy of a large sample of WLQs could provide measure¬ 
ments of these emission-line properties in the rest-frame op¬ 
tical. If they indeed have very high Eddington ratios, cluster¬ 
ing is expected toward high Fe ll/H/3, low [O III] REW, and 
low H/3 FWHM, relative to quasars with comparable lumi¬ 
nosities; there are already hints of such clusteri ng based on a 
small sample of five WLQs (see S ection 5.1 ofjPlotkin_et^ 
[20]^. Moreover, the Lliang et ^ (120141) and ISadowski et al. 
(120141) simulations of super-Eddington accretion disks predict 
a radiation-driven outflow along the rotation axis. Signatures 
of this outflow may be sought in our X-ray normal WLQs as 
high-ionization absorption lines/edges at X-ray and UV wave¬ 
lengths (e.g., O VI absorption, as discussed in Section 4.4 of 
Wll). 

To quantify the possible correlations between Aaox and 
Fe II REW or A(g - /) (Figure fTTl i and probe the underly¬ 
ing physics, deeper X-ray observations of the undetected ob¬ 
jects are required. Given the stacked X-ray flux level, most 
of the undetected sources should be detectable with factors of 
R 10 increase in the exposure times. Such observations will 
also help test whether the Aaox distribution is bimod al, a s 
might be expected in the shielding-gas scenario (Section lOT i. 
Moreover, long-term X-ray monitoring observations of the 
PHL 1811 analogs and WLQs will be useful, as they can con¬ 
strain the frequency and duration of the X-ray state transitions 
as in PHL 1092, and thus provide insights into the geometri¬ 
cally thick disk scenario. 

Finally, it remains somewhat perplexing that at least many 
of our PHL 1811 analogs appear to be X-ray absorbed, while 
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PHL 1811 itself appears intrinsically X-ray weak; Occam’s 
razor would initially favor a single explanation for the X-ray 
weakness of all these objects selected to have similar UV 
emission-line properties. It is worth noting that our PHL 1811 
analogs are generally being observed in the Ri 1.5-24 keV 
rest-frame band, while the X-ray properties of PHL 1811 it¬ 
self have only been effectively probed up to r; 8 keV. Perhaps 
PHL 1811 itself has a highly absorbed component that has 
yet to be recognized. Although PHL 1811 has not been de¬ 
tected in the 5w!/f-BAT all-sky survey in the 14-24 keV band 
(M. Koss 2014, private communication), a NuSTAR observa¬ 
tion could probe more sensitively for a highly absorbed com¬ 
ponent. 
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Table 1 

New Chandra Observations and X-ray Photometric Properties 


Object Name 

Redshift 

Observation 

Observation 

Exposure 

Soft Band 

Hard Band 

Band 

Teff 

Comment 

(J2000) 


ID 

Start Date 

Time (ks) 

(0.5-2 keV) 

(2-8 keV) 

Ratio 



(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 


PHL 1811 Analogs 


014733.58 + 000323.2 

2.031 

14949 

2013 Nov 15 

6.2 

<2.4 

3.2+” 

> 1.19 

<0.6 


082508.75 + 115536.3 

1.998 

14951 

2013 Jun 10 

5.1 

<2.4 

<2.5 




090809.13+444138.8 

1.742 

14950 

2013 Sep 23 

6.8 

<2.4 

<2.5 




094808.39+161414.1 

1.805 

14947 

2013 Jun 10 

3.7 

2 

I ^ 

<2.5 

< 1.42 

>0.4 


113342.67 + 114206.2 

2.052 

14952 

2014 Feb 08 

9.3 

< 2.4 

2 

^•^-1.4 

>0.70 

< 1.0 


143525.31+400112.2 
152156.48 + 520238.5 

2.267 

2.238 

14954 

15334 

2013 Sep 04 
2013 Oct 22 

6.9 

37.1 

43 5+’''* 

2 

47.9+*'5 

0.29:(“'« 

1 iq+I].28 
^•■^^-0.24 

1 52+1-3 

U.D-0 2 

Wll, archival 4 ks obs. 

153412.68 + 503405.3 

2.122 

14956 

2013 Dec 31 

6.1 

<4.1 

<4.2 




153714.26 + 271611.6 

2.457 

14955 

2013 Dec 22 

6.0 

44.4-« 

22 0+^'^ 

-0.13 

^•■^-0.3 


153913.47 + 395423.4 

1.934 

14948 

2013 Dec 13 

5.3 

<2.4 

<2.5 




222256.11-094636.2 

2.913 

14953 

2013 Sep 27 

9.4 


<4.2 

<2.87 

>-0.2 


Weak-line Quasars 

014333.65-391700.1 
082722.73 + 032755.9 

1.800 

2.031 

15358 

15342 

2012 Nov 21 

2013 Mar 23 

3.0 

2.5 

15 6+^'^ 
<2.4 

y y+4.2 
<2^1 

0.49+“:3^ 

1 4+0.5 
^ • -0.4 

HE 0141-3932 

084424.24+124546.5 

2.492 

15337 

2013 Apr 01 

1.5 

51.2+8-^ 

ll.O+t'l 

0.21+“'*“ 

2-+n1 


090843.25 + 285229.8 

0.933 

15352 

2013 Feb 25 

2.0 

40.6^1 

'- 2.8 

0 19’8:?I 
^•^^-0.08 

7 7 +[]:^ 

^•^-0.4 


100517.54 + 331202.8 

1.802 

15351 

2013 Mar 01 

3.0 

<4.1 

<2.5 




113413.48 + 001042.0 

1.487 

15343 

2013 Feb 19 

2.5 

<4.1 

<2.5 




115637.02+184856.5 

1.993 

15341 

2013 Mar 20 

2.0 

43 

11 0+^-^ 

A 75 + 0.12 

^•^-^-0.09 

2.o!j>:^ 


130312.89 + 321911.4 

0.638 

15350 

2012 Nov 19 

2.0 

<2.4 

<2.5 




131059.77 + 560140.2 

1.285 

15344 

2013 Apr 10 

2.0 


-’•J_2.4 

^•^^-0.14 

1.8!“:“ 


132130.21+481719.1 

1.409 

15347 

2013 Feb 19 

2.5 

<2.4 

<2.5 



132809.59 + 545452.7 
133222.62 + 034739.9“ 

2.116 

1.447 

15338 

15340 

2013 Jul 15 
2013 Feb 13 

2.0 

2.0 

6.2+” 

2 

7 7+2.9 
6:6’|| 

0 35+“'53 

^ ■ r o 1 

1 7+*'3 
+ ■'-0,8, 
-o.3!‘:> 

Lensed, BAL 

134601.28 + 585820.2 

1.664 

15336 

2013 Apr 28 

1.5 

<2.4 

<2.5 




140710.26 + 241853.6 

1.668 

15345 

2012 Nov 19 

2.5 

<2.4 

<2.5 




141141.96+140233.9 

1.753 

15353 

2012 Dec 16 

3.4 

48.9i” 

13.2!5:“ 

-0.09 



141730.92 + 073320.7 
142943.64 + 385932.2 

1.710 

0.928 

15349 

15335 

2012 Dec 05 
2012 Dec 08 

2.5 

1.5 

21.9!|'S 

14.6+5'“ 

6 2+3^ 

u.z 25 

-^•“^-2.4 

4.4t” 

^•^^-0.12 

Q 30-^^•26 
-0.16 

2.0!“'“ 

I'Ste 


144741.76-020339.1 

1.431 

15355 

2013 Jan 13 

2.0 

<4.3 

<0.70 

> 1.0 


162933.60 + 253200.6 

1.339 

15356 

2012 Dec 07 

2.9 

is.vii 

1 /-0^,4 

0.94!“:« 

0.8!“'^ 

2 1+0.4 
^■+-0.4 


164302.03+441422.1 

1.650 

15348 

2013 May 30 

3.0 

40.6;3« 

8.8^:3 



172858.16 + 603512.7 

1.807 

15346 

2013 Jan 05 

2.5 

^'2-^-4,0 

<4.3 

<0.29 

> 1.8 


212416.05-074129.8 

1.402 

15339 

2012 Nov 20 

1.5 

13.5+1'2 

<4.3 

<0.35 

> 1.7 


215454.35-305654.3 

0.494 

15354 

2012 Dec 13 

2.5 

<2.4 

<2.5 



2QZJ215454.3-305654 

Note. — Cols. (1) and (2): object 

name in the J2000 equatorial coordinate format and redshift. Cols. (3)-(5): 

Chandra observation ID, observation start date. 


and background-flare cleaned effective exposure time in the 0.5—8 keV band. Cols. (6)-(7): Aperture-corrected source counts in the soft (0.5-2 keV) and hard 
(2-8 keV) bands. An upper limit at a 90% confidence level is given if the source is not detected. Col. (8): ratio between the soft-band and hard-band counts. An 
entry of “...” indicates that the source is undetected in both bands. Col. (9): 0.5-8 keV effective power-law photon index, derived from the band ratio assuming a 
power-law spectrum modified with Galactic absorption. An entry of “...” indicates that it cannot be constrained. Col. (10): comments on special objects. 

^ J1332+0347 is included in this table but not in the figures or our analyses (see Section[22)- 
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Table 2 

Quasar UV Emission-Line Measurements from the SDSS Spectra 


Object Name 

MJD 

C IV Blueshift 

C IV FWHM 

REW 

REW 

REW 

REW 

REW 

REW 

(J2000) 


(km s“*) 

(km s“*) 

C IV 

Si IV 

A1900 A 

Fe II 

Fe III 

Mg II 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(V) 

(8) 

(9) 

(10) 

PHL 1811 Analogs 

014733.58+000323.2 

51793 

-6170 ±440 

7380 ±560 

4.5 ±0.3 

<3.0 

3.8±0.3 

21±2 

6.2±0.3 

12.2±0.6 

082508.75 + 115536.3 

54149 

-4600 ±340 

8100 ±640 

5.8±0.3 

4.0±0.3 

8.2 ±0.4 

26±2 

5.4±0.4 

13.1 ±0.8 

090809.13+444138.8 

52312 

-4250 ±640 

5650 ±460 

4.1 ±0.5 


7.9 ±0.5 

32±3 

4.8±0.8 

17.1±1.1 

094808.39 + 161414.1 

54095 

-3700 ±380 

8290 ±660 

5.9±0.3 


3.7 ±0.3 

28±2 

4.7 ±0.3 

16.5 ±0.5 

113342.67 + 114206.2 

53055 

-2460 ± 1740 

2900 ± 240 

2.6±0.8 

I.4±1.2 

4.9 ±0.9 

27±2 

3.4±0.7 

10.4±1.9 

143525.31+400112.2 

52797 

-2050 ±440 

6160 ±500 

7.4±0.6 

3.2 ±0.5 

8.3 ±0.8 

13±1 

1.9±0.5 


152156.48 + 520238.5 

52376 

-9300 ±610 

11700 ±800 

9.1 ±0.6 

2.7 ±0.3 

8.2 ±0.6 

21±2 

7.1 ±0.6 

14.9±0.6 

153412.68 + 503405.3 

52401 

-970 ± 320 

7610 ±580 

8.7 ±0.5 

3.1±0.3 

15.8±0.5 

40±4 

7.8 ±0.4 

25.6±1.0 

I537I4.26 + 27I611.6 

54180 

-1310±220 

6560 ±520 

6.0 ±0.3 

3.9 ±0.2 

8.0±0.3 

> 11 (13±2) 

1.9±0.3 


153913.47 + 395423.4 

53171 

-3930 ±600 

10570 ±740 

6.9 ±0.5 

2.2 ±0.4 

11.2±0.6 

30±3 

5.4±0.4 

19.9±0.8 

222256.11-094636.2 

52206 

-2570 ± 1300 

5030 ±440 

3.7±0.5 

5.2 ±0.9 

8.7 ±0.7 


< 3.0 


Weak-line Quasars 

014333.65-391700.1 

082722.73+032755.9 

52642 

-1460 ±520 

2560 ±200 

4.0±0.8 

< 1.1 

2.5±1.1 

53±5 

6.0±0.8 

22.2 ±2.1 

084424.24+124546.5 

53801 

-1450 ±500 

5870 ±480 

5.5±0.5 

2.1 ±0.4 

3.1±0.4 

> 16 (26 ±2) 

3.9±0.5 


090843.25 + 285229.8 

53330 






19±1 

< 1.6 

4.9±0.8 

100517.54 + 331202.8 

53378 

-2420 ± 1040 

4020 ± 340 

5.8± 1.0 


3.4±1.0 

24 ±2 

2.6±1.3 

12.5±2.1 

113413.48+001042.0 

51630,51658 





5.1±0.5 

38±3 

3.2±0.5 

12.9±0.6 

115637.02 + 184856.5 

54180 

-1790 ±900 

3180 ±260 

1.2±0.5 

1.5 ±0.4 

5.0±0.8 

22±2 

1.8±0.6 

8.5±1.5 

130312.89 + 321911.4 

53819 






> 10 


7.3±0.9 

131059.77 + 560140.2 

52791 





4.5 ±0.8 

12±1 

1.2±0.5 

7.3±1.0 

132130.21+481719.1 

52759 





5.6±1.0 

27 ±2 

3.0±0.9 

9.3±1.0 

132809.59 + 545452.7 

52724 

-2560 ±340 

3970 ± 300 

4.4 ±0.4 

3.5 ±0.3 

3.7 ±0.5 

20±2 

3.5±0.5 

18.3±2.2 

133222.62+034739.9"' 

52374 





4.4±1.1 

8±1 

< 1.5 

9.3±0.9 

134601.28 + 585820.2 

52425,52466 

-4570 ±760 

5580 ±440 

2.9±0.3 


4.8 ±0.3 

14±1 

1.7±0.3 

12.1 ±0.5 

140710.26 + 241853.6 

53770 

-780 ±600 

1170 ±140 

1.5±0.6 


5.4±0.8 

26±2 

4.1 ±0.9 

10.0±1.2 

141141.96+140233.9 

53442 

-1850 ±740 

2420 ±180 

2.4±0.9 


6.6 ±0.9 

27 ±2 

5.5±1.2 

6.0 ±1.6 

141730.92+073320.7 

53499 

-5470 ±760 

3240 ± 260 

2.0±0.7 


4.2 ±0.8 

25±2 

2.8±0.8 

11.2±1.2 

142943.64 + 385932.2 

52797 






22±2 

2.3±0.5 

13.7±0.8 

144741.76-020339.1 

52411 





6.7 ±0.7 

7±1 

2.8±0.7 

12.4±1.1 

162933.60 + 253200.6 

53226 





9.5 ±1.0 

5±1 

3.2±0.9 

10.7±1.1 

164302.03+441422.1 

52051 



< 2.0 


<2.0 

8±1 

9.8±1.6 

3.9±1.1 

172858.16 + 603512.7 

51792 

-2630 ±720 

3380 ±340 

3.1 ±0.7 

2.6±1.1 

6.2 ±0.8 

11±1 

< 1.6 

5.9±1.6 

212416.05-074129.8 

215454.35-305654.3 

52178,52200 





2.6±0.5 

25±2 

1.4±0.4 

12.7±0.5 


Note. — Col. (1): object name. Col. (2): modified Julian date of the SDSS observation; an entry of “...” indicates that no SDSS spectiirm is available. In 
a few cases, two spectra are available, and the average spectrum was used for the measurements. Cols. (3) and (4): blueshift and full width at half maximum 
(FWHM) of the Civ line. Cols. (5H10): REWs (in units of A) of the C IV A1549, SilV A1397, A1900 complex, Fe II (2250-2650 A), Fe III UV48 A2080, 
and Mg II A2799 emission features. The Si IV A1397 line is a blend of Si IV and O IV], and the A1900 emission is domi nated by C III] A1909 wit h additional 
contributions from [Ne III] A1814, Si II A1816, A1 III A1857, Si III] A1892, and several Fe III multiplets (see Table 2 of IVanden Berk et alJ2001l for details). 
The UV Fe II REW was measured between 2250A and 2650A; for three objects that do not have full spectral coverage of this range, a lower limit was derived 
(along with an estimated value in parentheses if the covered fraction is > 60%). For the other properties, the measurements were performed the same manner 
as in Wll. The uncertainties have been multiplied by the factors given in the notes of W11 Table 3, and the upper limits ai‘e at a 3cr confidence level. 

^ J1332-I-0347 is included in this table but not in the figures or our analyses (see Section lT^ . 














Table 3 

X-ray and Optical Properties 


K) 


Object Name 

Mi 

A'H.Gal 

Count Rate 

Px h keV 

logix / 2500 A 

log 7-2500 A 

Qox 

Aq!ox(o') 

/weak 

Mg-i) 

R 

(J2000) 



(0.5-2 keV) 

(0.5-2 keV) 

(2-10 keV) 






(1) 

(2) 

(3) 

(4) 

(5) (6) 

(7) (8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 


PHL 1811 Analogs 


014733.58 + 000323.2 

-27.91 

2.90 

<0.39 

<0.16 

<0.68 

< 43.77 

3.68 

31.58 

< 

-2.20 

< 

-0.49(3.32) 

> 18.40 

0.55 

< 

1.0 

082508.75 + 115536.3 

-28.31 

3.67 

<0.47 

< 0.19 

< 0.81 

<43.83 

4.62 

31.66 

< 

-2.21 

< 

-0.48(3.30) 

> 17.95 

-0.08 

< 

1.0 

090809.13+444138.8 

-27.26 

1.66 

<0.35 

< 0.14 

< 0.56 

<43.56 

2.09 

31.21 

< 

-2.14 

< 

-0.47(3.25) 

> 17.28 

0.00 

< 

1.9 

094808.39+161414.1 

-28.43 

3.34 

0 55+0.74 

0.22 

0.90 

43.80 

6.59 

31.74 


-2.25 


-0.51(3.51) 

21.67 

0.20 

< 

0.6 

113342.67 + 114206.2 

-27.52 

3.32 

< 0.26 

<0.10 

< 0.45 

< 43.60 

1.99 

31.32 

< 

-2.17 

< 

-0.49(3.34) 

> 18.61 

0.39 

< 

2.3 

143525.31+400112.2 

-27.84 

1.00 

I 057^-^’^ 
■^-0.39 

0.42 

1.89 

44.32 

3.32 

31.62 


-2.01 


-0.29(2.00) 

5.75 

0.11 

< 

1.4 

152156.48 + 520238.5 

-30.22 

1.58 

1 17+<’'21 

+ ■+Ig 

0.46 

0.98 

44.51 

19.86 

32.39 


-2.42 


-0.59(4.51) 

34.73 

0.24 

< 

0.2 

153412.68 + 503405.3 

-28.70 

1.56 

< 0.67 

< 0.27 

< 1.20 

< 44.05 

3.69 

31.61 

< 

-2.11 

< 

-0.39(2.64) 

> 10.13 

0.49 

< 

1.2 

153714.26 + 271611.6 

-28.67 

3.06 

7 45+1.32 

2.92 

10.23 

45.31 

5.75 

31.92 


-1.82 


-0.06(0.45) 

1.43 

0.03 

2 

;.6 

153913.47 + 395423.4 

-27.99 

1.70 

< 0.45 

<0.18 

<0.76 

< 43.78 

3.86 

31.56 

< 

-2.19 

< 

-0.48(3.26) 

> 17.39 

0.09 

< 

1.2 

222256.11-094636.2 

-28.55 

4.34 

0 224^-^^ 

0.09 

0.47 

43.88 

4.33 

31.92 


-2.29 


-0.52(4.01) 

23.29 

0.23 

< 

1.3 


Weak-line Quasars 


014333.65-391700.1 

-29.22 

1.71 

5.28+{'’^ 

2.07 

6.75 

44.88 

15.28 

32.10 


-2.06 

-0.27(2.03) 

4.93 


< 0.5 

082722.73 + 032755.9 

-27.27 

3.73 

< 0.97 

< 0.39 

< 1.68 

<44.16 

1.68 

31.24 

< 

-1.92 

<-0.25(1.72) 

>4.50 

0.31 

<2.6 

084424.24+124546.5 
090843.25 + 285229.8 

-28.17 

-25.10 

3.81 

2.46 

33.58!|'g 

20.28+|g 

13.53 

8.30 

74.79 

23.71 

45.90 

44.55 

3.81 

1.48 

31.75 

30.54 


-1.42 

-1.46 

0.32(2.17) 

0.11(0.57) 

0.15 

0.51 

0.37 

-0.03 

8.0 

<2.0 

100517.54 + 331202.8 

-26.70 

1.45 

< 1.38 

<0.55 

<2.22 

<44.19 

1.38 

31.06 

< 

-1.84 

<-0.20(1.35) 

>3.26 

0.07 

<2.8 

113413.48 + 001042.0 

-26.48 

2.63 

< 1.64 

<0.66 

<2.39 

< 44.07 

1.57 

30.96 

< 

-1.85 

<-0.22(1.11) 

>3.75 

0.17 

<2.4 

115637.02+184856.5 

-27.30 

2.83 

21.91C3“ 

8.78 

38.74 

45.49 

2.04 

31.31 


-1.43 

0.25(1.71) 

0.22 

0.05 

7.4 

130312.89 + 321911.4 

-24.20 

1.17 

< 1.20 

<0.48 

< 1.20 

< 43.03 

1.45 

30.19 

< 

-1.95 

<-0.43(2.17) 

> 13.09 

0.21 

< 1.9 

131059.77 + 560140.2 

-26.09 

1.51 

9.36+2-]’ 

3.73 

12.42 

44.69 

1.64 

30.86 


-1.58 

0.03(0.17) 

0.82 

-0.10 

<2.3 

132130.21+481719.1 

-26.20 

1.13 

< 0.97 

<0.39 

< 1.36 

< 43.78 

1.38 

30.86 

< 

-1.92 

<-0.31(1.54) 

>6.26 

0.12 

<2.6 

132809.59 + 545452.7 

-27.72 

1.16 

Q 1 7+1.86 

1.23 

4.89 

44.75 

2.70 

31.47 


-1.82 

-0.12(0.81) 

2.03 

-0.09 

< 1.7 

133222.62 + 034739.9“ 

-26.70 

1.98 

1 57 

0.41 

0.63 

44.47 

1.75 

30.98 


-2.09 

-0.46(2.31) 

15.58 

0.34 

< 2.2 

134601.28 + 585820.2 

-27.38 

1.63 

< 1.57 

<0.63 

<2.43 

<44.17 

3.22 

31.36 

< 

-1.97 

<-0.28(1.92) 

>5.37 

0.11 

< 1.3 

140710.26 + 241853.6 

-26.70 

1.65 

<0.97 

<0.39 

< 1.50 

< 43.96 

1.50 

31.03 

< 

-1.92 

<-0.28(1.41) 

>5.36 

0.14 

<2.8 

141141.96+140233.9 

-26.56 

1.43 

14.42+|:« 

5.75 

22.89 

45.18 

1.13 

30.95 


-1.42 

0.21(1.06) 

0.28 

-0.03 

<3.6 

141730.92 + 073320.7 

-26.64 

2.12 

8.81!2|’ 

3.53 

14.14 

44.93 

1.64 

31.09 


-1.56 

0.09(0.60) 

0.59 

-0.01 

<3.0 

142943.64 + 385932.2 

-26.01 

0.95 

9 56++“® 
’■■“'’-2.53 

3.81 

11.02 

44.37 

3.57 

30.91 


-1.73 

-0.11(0.55) 

1.91 

-0.05 

<0.9 

144741.76-020339.1 

-26.00 

4.53 

Q 1 1 +T.86 

■^•-^^-1.24 

1.25 

4.46 

44.31 

1.69 

30.96 


-1.76 

-0.13(0.64) 

2.15 

-0.36 

<2.3 

162933.60 + 253200.6 

-25.35 

3.58 

6.39!j-| 

2.51 

6.44 

44.87 

1.05 

30.70 


-1.62 

-0.02(0.12) 

1.16 

0.01 

<3.4 

164302.03+441422.1 

-26.56 

1.52 

13.73+2;]’ 

5.52 

22.46 

45.06 

1.20 

30.93 


-1.43 

0.19(0.98) 

0.31 

0.13 

<3.1 

172858.16 + 603512.7 

-26.88 

3.32 


2.52 

10.21 

44.86 

1.86 

31.19 


-1.64 

0.03(0.18) 

0.85 

0.05 

<2.4 

212416.05-074129.8 

-26.64 

5.58 

8.89!+2‘ 

3.60 

12.68 

44.74 

2.56 

31.12 


-1.65 

-0.00(0.00) 

1.00 

-0.14 

< 1.5 

215454.35-305654.3 

-23.19 

1.84 

< 0.97 

<0.39 

< 0.90 

< 42.68 

0.61 

29.59 

< 

-1.85 

< -0.42(2.53) 

> 12.25 


<8.3 


Note. — Col. (1): object name. Col. (2): absolute /-band magnitudes. Col. (3): Galactic neutral hydrogen column density iDickev & Lockman|199QI) . Col. (4); observed 0.5-2 keV Chandra count rate in 
units of 10~^ s~^ Col. (5): Galactic absorption-corrected observed-frame 0.5-2 keV flux in units of 10”^“^ erg cm~^ s~^ Col. ( 6 ): rest-frame 2 keV flux density in units of 10~^^ erg cm~^ s~^ Col. (7): 

logarithm of the rest-frame 2-10 keV luminosity in units of erg s~^, derived from the observed 0.5-2 keV flux. Col. (8); flux density at rest-frame 2500 A in units of 10~^^ erg cm~^ s~^ . C ol. (9): logarithm 

of the rest-frame 2500 A luminosity density in units of erg s~^ Hz~^ Col. (10): measured chqx pai'ameter. Col. (11): difference between th e measured chqx and the expected aox from the lJust et alll20d% 
q:ox-^2500 a elation. The statistical signiflcance of this difference, measured in units of the qqx rms scatter in Table 5 of ISteffen et all [ 200^, is given in the p arenthesis. Col. (12): factor of X-ray weakness in 
accordance with Achqx- Col. (13): relative SDSS g-i color. Col. (14): radio-loudness parameter, defined as /? = /s nw// i_ (e.g.. lKellermann et all 19891) . where fs ghz and /4400 A densities at 

rest-frame 5 GHz and 4400 A, respectively. 

^ J1332+0347 is included in this table but not in the figures or our analyses (see Section fZ^ . 
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Table 4 

Stacked X-ray Properties for the X-ray Weak Subsamples 


Subsample ] 

Number of 
Sources 

Mean 

Redshift 

Total Stacked 
Exposure (ks) 

Soft-Band 

Counts 

Hard-Band 

Counts 

Teff 

aox 

AqoxCo-) 

P1811A‘' 

15 

2.16 

106.6 

26.5!|'4 

16.l1« 

1 1540.37 

- 2.33 

-0.61(4.16) 

P1811Aundet‘> 

7 

2.08 

46.8 

2.9+f“ 

18 

Q 79+1.63 
\j.iy Q g3 

-2.57 

-0.87(5.94) 

PlSllASBundet" 

9 

2.07 

62.3 

2.8":* 

7 9+4.^ 

-0 

-2.74 

-1.04(7.09) 

WLQ‘' 

15 

1.53 

61.7 

31-811 

i2-2l:§ 

1 ^c+0.40 
^•^°-0.35 

-2.15 

-0.47(3.22) 

WLQ undet*^ 

10 

1.54 

31.4 

i-8l:I 

< 2.5 

> 0.44 

-2.44 

-0.78(5.35) 

C IV WLQf 

7 

1.89 

36.8 

5.0+|j 

< 3.8 

> 1.1 

-2.30 

-0.62(4.27) 

C IV WLQ undet^ 

6 

1.89 

21.9 

<3.9 

< 2.5 


< -2.26 

< -0.59(4.02) 

Allh 

30 

1.85 

168.3 

58-4l:§ 

28.3+’;“ 

1 7J7+0.25 

-0.23 

-2.24 

-0.54(3.68) 


^ X-ray weak PHL 1811 analogs with Chandra observations except J1521-1-5202 (nine in the Chandra Cycle 14 sample and six in 

Wll). 

^ Undetected PHL 1811 analogs with Chandra observations (four in the Chandra Cycle 14 sample and three in Wll). 

Soft-band undetected PHL 1811 analogs with Chandra observations (six in the Chandra Cycle 14 sample and three in W11). The are 
two objects that are not detected in the soft band but are detected in the hard band. This stacked Aqqx value was used in Figure[^. 

^ X-ray weak WLQs with Chandra observations (nine in the Chandra Cycle 14 sample and six in W12). 

^ Undetected WLQs with Chandra observations (eight in the Chandra Cycle 14 sample and two in W12). This stacked Aaox value 
was used in Figure[6j). 

^ X-ray weak WLQs with Chandra observations in the C IV subsample (four in the Chandra Cycle 14 sample and three in W12). 

® Undetected WLQs with Chandra observations in the C IV subsample (four in the Chandra Cycle 14 sample and two in W12). 

^ All X-ray weak PHL 1811 analogs and WLQs with Chandra observations except J1521-I-5202. 


Table 5 

Bolometric Luminosity, SMBH Mass, and Eddington Ratio 
Estimates 


Object Name 
(J2000) 

logLBol 
(erg s“*) 

logMBH 

(M©) 

+Bol/+Edd 

Chandra Cycle 14 PHL 1811 Analogs 

014733.58+000323.2 

46.8 



082508.75 + 115536.3 

47.5 



090809.13 + 444138.8 

46.8 

9.8 

0.07 

094808.39 + 161414.1 

47.5 

9.9 

0.38 

113342.67 + 114206.2 

46.8 



143525.31+400112.2 

46.8 



152156.48+520238.5 

48.2 

9.8* 

2.09 

153412.68 + 503405.3 

46.8 



153714.26 + 271611.6 

46.8 



153913.47 + 395423.4 

47.3 




Note. — The SMBH masses are the M g Il-based single¬ 
epoch virial masses from IShen et alJd201 II) . except for those 
objects marked with a which are th e H/3-based virial 
masses from Wll or lPlotkin et all I2015D . We caution that 
the SMBH mass and bolometric lum inosi ty estim ates have 
substantial uncertainties (see Sections [5.3l and [4!^ . and thus 
so do the Eddington ratios. 

Only a portion of this table is shown here to demonstrate its 
form and content. A machine-readable version of the full 
table is available. 
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Table 6 

Peto-Prentice Test Results for the Spectral Properties of X-ray Weak and X-ray Normal 

Samples 


Spectral 

Property 


Full Sample 



C IV Subsample 

•^weak 

^normal 

<7 


•^weak 

^normal 

(7 

^llull 

CivREW 

25 

11 

2.3 

0.02 

25 

11 

2.3 

0.02 

C IV blueshift 

23 

9 

1.7 

0.09 

23 

9 

1.7 

0.09 

C IV FWHM 

23 

9 

2.2 

0.03 

23 

9 

2.2 

0.03 

A1900 ArEW 

29 

16 

1.7 

0.10 

24 

11 

1.7 

0.10 

Fe II REW 

30 

18 

3.8 

1 X lO^'* 

24 

11 

2.9 

4 X 10-3 

Fe III REW 

29 

18 

2.5 

0.01 

24 

11 

1.2 

0.21 

Mg II REW 

21 

16 

2.4 

0.02 

15 

9 

2.7 

8 X 10-3 

Mg-i) 

31 

18 

4.6 

4 X 10-'’ 

25 

11 

3.7 

2x 10^ 


Note. — In cases where there are no censore d data, the Peto-Prentice test reduces to 
Gehan’s Wilcoxon test (e.g.. lLavallev et ^I1992ll . The tests were performed for both the 
full sample and the C IV subsample. We list for each test the number of X-ray weak objects 
(A/weak)^ number of X-ray normal objects (Abnormal)’ t^st statistic (cr), and probability of the 
data being drawn from the same parent population (/null)- 










